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ABSTRACT 


A Boolean function may be expressed, using Exclusive 
OR, AND, and NOT gates, in an unique canonical form, called 
Reed-Muller expansion, which is similar to minterm expansion 
of a Boolean function. The expansion with Exclusive-QR gates 
possesses inherent advantages over conventional realisation 
of Boolean functions. Minimisation of Exclusive-OR gates, 
to suit optimal design conditions, consists of obtaining what 
is known as the polarity function of the coefficients in Reed- 
Muller expansions; polarity function give^. logical conditions 
of the input variables. Transformations of minterm coeffi- 
cients to the coefficients in Reed-Muller expansions^being 
linear, are represented by matrices, called Reed-Muller 
matrices. Reed-Muller matrices are factorised using direct - 
product £•■ Good’s techniques, which facilitates handling of 
these matrices in the computer. Algorithmic methods for 
finding the eigen-structure of Reed-Muller matrices, using 
the concepts of the null-space, generalised inverses of 
matrices, have been evolved and : . utilized for the genera- 
tion of coefficients in Reed-Muller expansions from which 
one can obtain polarity functions for the minimum 
expansion. 
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CHAPTER 1 


INTRODUCTION 


1.1 SCOPE OF WORK: 

An attempt has been made in this thesis to study 
the eigenstructure of Reed-Muller matrices and utilize 
their structure for generation of Reed-Muller expansion 
coefficients of Boolean functions. An alternative procedure 
based on factorization of Reed-Muller matrices, using Good's 
techniques, has also been given for the 
generation of Reed-Muller expansion coefficients. 

Any Boolean function may be expressed in a sum(s)- 
of-product canonical form. Employing only Exclusive-OR 
(Ex-OR) gate and AND gates and by algebraic substitution and 
simplification this may be uniquely represented by a similar 
canonical form [l],[2]. As an example, a 3 variable function 
of (x 3 ,x 2 ,x 1 ) is considered. The m interm expansion , 
using AND gates, OR gates and NOT gates of this function is 


f ( x 3 ,x 2 ,X l) -C o X 3 X 2 X 1 + C 1 X 3 X 2 X 1 + c 2 X 3 X 2 X 1 + c 3 X 3 X 2 X 1 


+ c 4 x 3 x 2 x i +c 5 x 3 x 2 x 1 + c 6 x 3 x 2Xl + c ? x 3 x 2 x 2 


( 1 . 1 ) 



2 


where c. fo,l? depending on the function value of the related 
input variables. Every OR gate in the above expression may 
be replaced by an Ex-OR gate as only one minterm becomes 
1 for every combination of values x 3 ,x 2 ,x^ and others 0* The 
equivalent expression, therefore, may be written as 

|-( x 3» x 2» x i) = c 0 *3 X 2 X 1 © c i X 3 X 2 X 1 © c 2 X 3 X 2 X 1 C 3 X 3 X 2 X 1 

t) 

© c 4 X 3 X 2 X 1 © c 5 X 3 X 2 X 1 © c 7 X 3 X 2 X 1 

( 1 . 2 ) 

* 

and substituting the identity = 1 © *i> the expansion 
is obtained in the following form 

i c (x 3 ,x 2 , x i) = c o (l@ x 3 )U©x 2 )(l© x 1 )© c 1 (l©x 3 ) 

(1©x 2 ) x i© c 2 (l ® x 3 )x 2 (l©x 1 ) © c 3 (l©x 3 ) 

X 2 X 1© C 4 x 3 (1 © x 2 )(1 © x l ) © c 5 x 3 (1 © x 2 ) X 1 
@ c 6 X 3 X 2^ 1 © X 1^ < © c 7 x 3 x 2 X 1 

(1.3) 

which reduces to 

= c 0 @ (c o© c l )x l© (c o© C 2 ) x 2 
0( c o© c l© c 2© c 3 )x 2 x l© (c o® c 4 ) x 3 
©( c 0 © c l© c 4© c 5^ X 3 X 1©) 
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( c o © c 2® C 4@ c 6^ X 3 X 2©^ C o® C l© 
c 2© c 3© c 4 © c 5® c 6 © C 7^ X 3 X 2 X 1 
~ a 0 © a l x l © a 2 X 2 © a 3 x 2 X l © a 4 x 3© 

a 5 x 3 x l@ a 6 x 2© a 7 x 3 x 2 x l (l * 4) 


where 

a o * c o 

a l = c o© c l 

a 2 * c o & c 2 

a 3 = c 0 0 C 1 8 c 2 © c 3 

a 4 = C o © c 4 

a 5 = c 0 © c i €) c 4 © c 5 

a 6 = C o© C 2© C 4© c 6 

a 7 » c 0 © c i® c 2© c 3© c 4© c 5® c 6© c 7 


(1.5) 

and a. ^0,lw are constant coefficients that are determined 

1 L J 

by the value of the function for each combination of values 
Xg, x 2 » x i* As thereto complemented variables in (1.4), 
this is also termed as a polynomial expansion. This expansion 
is unique because the expansion with Ex-OR and minterm is 
unique and the conversion to the expansion with a.s as the 
coefficients is also unique. 
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Thus any m variable Boolean function may be expressed 
as a polynomial expansion. These are also known as positive - 
canonic or complement free Reed-Muller (CFRM) expansions of 
2 m terms [3], [4], In general, the expansion may be written as 


f ( x l> x 2* “ a o ® a l x l © a 2 x 2® * *• a 2 m -l x l x 2 **’ x m 


( 1 . 6 ) 

By complementing any variable a non-positive canonic expansion 
also called as the generalised Reed-Muller (GRM) expansion 
may be obtained, and is written as 


f (x 


l' x 2 ,x 3* ** 


* m ) * a o@> a i k i& 3 2*2® 


a 2 m -l X l x 2 X 3* * * x m 


(1,7) 


where each x^ input variable may appear in positive or negative 
form but not in both forms. 

Hence, an m variable Boolean function may possess 2 m 
different positive and non-positive canonic RM expansions, 
depending on the polarity of the input variables, which may 
be called as the polarity function of the variables of the 
function. Using the following identities it is' possible to - 
obtain any non-positive canonic expansion from either the 
positive canonic or any other non-positive canonic RM 
expansions. 
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x i = i (i) x i and x i = i® x i 

1@1=0 x i @x i =0 (1.8) 


Suppose the positive canonic expansions in 3 variables is 
given as 


f(x 3 ,x 2 ,x 1 ) a a o ©a 1 x 1 ©a 2 x 2 ® a 3 x 2 x i ® a 4 x 3 

CD a 5 x 3 x j a 6 x 3 X 2^^ a 7 x 3 x 2 X l ^ 1»9) 


then the expansion with polarity function 110 may be ootained 
by the substitution of x^ = 10Xp and therefore obtain 

f(x 3 ,x 2 »x 1 ) = (a 0 ©a 1 )@ 3^0(32 033) x 2 ©a 3 * 2 * 1 

a> 

CD ( a 4© a 5) x 3^© a 5 x 3 x i©( a 6(3) a 7^ X 3 X 2^ 
a 7 X 3 X 2 X 1 


(1.10) 


coefficients 

Where a^s axe the £ of the positive canonic RM 

expansions. Allowing the coefficients of this expansion 
to be de-noted as , (the decimal superscript devoting the 
decimal value of the binary, which represents the variables 
being complemented) these coefficients may be expressed in 
terms of the minterm. coefficients as 
a o l = %© a x = c x 

a l l = a i = c 0 © c i 
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*2 1 = a 2@ a 3 
a 3 X = a 3 
a 4 X * a 4© a 5 
■j 1 = a 5 
a 6^ = a 6© a 7 

Z3 


“ °1 ® c 3 

= c 0 © c i© c 2© c 3 

= c l© c 5 

= c 0 © c 1® c 4©°6 
= c 1© c 3© C 5© C 7 
= c 0 ® c l®' : 2®) C 3® c A 
©=5© c 6© c 7 

( 1 * 11 ) 


The constant coefficients a^s of the positive/non- 
positive canonic RM expansion are respectively called as 
positive/non-positive canonic RM expansion coefficients ox 
RM/GRM expansion coefficients also. The constant coefficients 

c i s, of the minterm expansion is accordingly, called as the 

or 

minterm expansion coefficient merely minterm coefficients. 


The classical problem in switching theory has been 
the realisation of a switching function with minimum number 
of logical circuit elements, Karnaugh map and Quine-TMcClusky* s 
algorithms have helped obtain a limited minimum conditions. 
Besides these algorithms, decomposition and other techniques 
have been employed to achieve economical multilevel expressions 
for switching functions [5], [6]. Geometric formats like 
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the mapping of the positive-canonic RM expansion coefficients 
are used to generate any Ex-OR realisation [7] and thus try 
and obtain a minimum Ex-OR design. Boolean realisations use 
AND/OR/NAND and NOR gates, whereas these same functions can 
also be implemented using Ex-OR/AND/NOT gates. 

The advantages of Ex-OR logic design are that more 
economical realisation in terms of the number of gates and 
the number of gate interconnections are possible; also, designs 
with Ex-OR gates facilitate easier testing in comparison 
with vertex networks since the change of any input to an Ex-OR 
gate will propogate a change to the output, unlike the 
the vertex Boolean realisation which require specific input 
pattern changes to sensitise a path to the output [8]. 

However, the testability advantage may even exceed the 
possible disadvantage of the cost factor in the Ex-OR 
realisation. 

The mapping of the minterm coefficients onto the 
RM expansion coefficients may be considered as the mapping 
of the vector space C to a vector space A over the binary 
field F. As this mapping may be written as a set of n equations 
in n variables viz. the a i s and where n = 2 111 , m being the 
number of variables in the Boolean- function, the mapping may 
be considered as a linear transformation T:C— ^A, which 
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satisfies the condition (g c x (+)h c ) T = g( c x T) (J)h( c^T) , 
fox all vectors c x and c in C and all scalars g,h in F, viz. 
0,1. If c x> c 2» #, * c n is any basis of the vector space C and 
a l ,a 2*** a n are an ^ vec ^ ors A, then there is one and only 
linear transformation T:C-^A with = a^.^c T = a n [ 9 ]. 
Allowing i^ to be the unit vectors, i^ = (1,0,... 0), ... 
i n = (0,0,... 1) of the vector space G, then each may be 
written as 

h T = a i - < r u' r i2 r m> 

1 2 T = a 2 = (*21 ,r 22» •••' r 2n ) 


i T = 
n 


n 


* (x. 


nl' n2 


> • • 


* r nn } 


( 1 . 12 ) 


and therefore there is just one linear transformation 

associated with (1.12) and this transformation is determined 

by the nxn matrix R = (r^). Therefore there is a one-to-one 

correspondence between the linear transformation T:C-^A and 

the nxn matrix R, with entries from the binary field F. 

Given T, the corresponding matrix R is the matrix with ith 

row, the row of co-ordinates i.Tj given R=(r. .), T is the 

3- X j 

unique linear transformation carrying each unit vector i^ 
of C into the ith row ( r jj/ r i2* # * r ^ n ) A. 
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Therefore, the processing of the minterm coefficients 
c^s may be expressed in the matrix form; by allowing the 
minterm coefficients and the GRhA expansion coefficients to be 
column vectors; this may be written as 

a = Rc (1.13) 

where R is a constant matrix, called as the Reed-Muller 
matrix, with matrix multiplication operation being modulo 2, 
as all operations are over GF(2). The example of the 3 
variable Boolean function may be written as 


a o 


. i 

0 

0 

0 

0 

0 

0 

0°1 


X 1 

0 i 

a l 


i 

i 

0 

0 

0 

0 

0 

0 


^ C 1 ! 

a 2 


i 

0 

1 

0 

0 

0 

0 

0 


C 2 

a 3 


' i 

1 

1 

1 

0 

0 

0 

0 


c 3 

a 4 


i 

0 

0 

0 

1 

0 

0 

0 

i 

• 

i 

C 4| 

a 5 


i 

1 

0 

0 

1 

1 

0 

0 : 

* * i 

1 

C 5 j 

a 6 

! 

• i 
« 

0 

1 

0 

1 

0 

1 

0 j 

i 

! 

! 

C M 

; a 7 


i 

1 

1 

1 

1 

1 

1 

ij 

< 


c o ® c i 

c o© c 2 

c o© c l© c 2© c 3 
c o© c 4 

C o© C l© C 4© c 5 

c o© c 2<5) C 4© c 6 

S 0 © c i© c 2© c 3© c 4© c 5© c 6© C 7 ' 
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The RM matrix for the positive canonic RM expansion 
may also be iteratively built up from [lO]. 



[ 1 ] 


(1.14) 


where m 1 = m-1, for the m variable Boolean function. 

In a similar fashion, it must be feasible to obtain 
the RM matrices for the non-positive canonic RM expansions. 

The scheme developed is not an iterative one as in ( 1.14) 
but the matrix is obtained as a direct-product of two kinds 
of matrices of order 2x2. These matrices are termed as core 
•matrices. The core matrix associated with a positive polarity 
variable is the 2x2 matrix and the one associated with 

the complemented variable is *] . RM matrices for the 

non— positive canonic RM expansions are thus obtained as a 
direct-product of a combination of these two kinds of core - 
matrices, combined as per the polarity function of the 
variables in the Boolean function. RM matricies may also be 
factorised using Good’s technique. 

Based on the manner in which the RM matrices may be 
factorised two classes of RM matrices may be had viz,, the 
RM matrix for the positive canonic RM expansion and for the 
non-positive canonic RM expansions. Salient features of the 
RM matrices observed are: 
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i) the RM matrix for the positive canonic RM expansion 
is involutary. 

ii) the roots of the characteristic equation of the RM 
matrix, for the positive canonic RM expansion, is 
in the GF(2) 

iii) the roots of the characteristic equation of the RM 
matrix is repetative 

iv) the characteristic equation of the RM matrix, for the 
non-positive canonic RM expansions, is a product of 
second degree irreducible polynomial, 

v) the roots of the characteristic equation” of the RM 
matrix, for the non-positive canonic RM expansions, 
are therefore to be found in the extension field 
also, 

vi) the number of kinds of roots being limited to the 
elements of the extension field, the roots of the 
characteristic equation of the RM matrix, for the 
non-positive canonic RM expansion, is also repetitive. 

The characteristic roots are found from the characteri 
stic equation of the matrices, and therefore, if characteri- 
stic roots are obtained as direct product of eigenvalues of 
the core-matrices, from which the RM matrices are constructed 
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then it is possible to construct the characteristic equation 
from these characteristic roots. The characteristic roots are 
known as eigen-values also. 

Every eigen-value has an eigen-vector corresponding 
to it. When the roots of the characteristic equations are 
repetitive, unique eigen-vectors for each repetitive root 
can not be found by the normal procedures. Two methods of 
obtaining the modal matrix, the matrix whose columns are the 
set of linearly independent (LIN) eigen-vectors of a matrix, 
have been outlined. In one, the concept of generalised 
inverses of matrices is used, and m the other the concept 
of null-space of a matrix is employed. The second method has 
been suitably adopted, for the RM matrix for the positive canonic 
RM expansion, and a simplified procedure has been evolved. 

Having obtained the eigen structure of the RM matrices, 
algorithms to generate the expansion coefficients from the 
given truth-table of a Boolean function, have been evolved. 

The algorithm for the positive -canonic RM expansion makes use 
of the eigen-vectors of the RM matrix and relates all 
possible minterm coefficients to its corresponding expansion 
coefficients. The set of eigen-vectors forms a subgroup of 
the group of all possible minterm coefficients and therefore, 
cosets may be formed and these minterm coefficients may 
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now be associated to the corresponding expansion coefficients 
through the eigen-vectors. The other algorithm makes use of 
the eigen-values for generation of the non-positive canonic 
RM expansion coefficients. Good’s technique for factori- 
sation of RM matrices may also be used for the generation of 
the RM expansion coefficients. 

A geometric mapping of the positive canonic RM 
expansion coefficients, called as 'b. maps’ , has been 

J 

developed to generate all non-positive canonic RM expansion 
coefficients [ 7 ]. However, the map operations become 
difficult when the number of variables m the Boolean function 
increases beyond 6. A computer method has been developed, 
to generate all non-positive canonic RM expansion coefficients 
from the positive canonic RM expansion coefficients [lO]. 

1.2 LINE OF APPROACH: 

The RM matrices were defined in the previous section ithe 
RM matrix are constructed as direct product of core-matrices. 
Another method of factorising the RM matrix based on 
Good’s technique is also considered. They are classified 
as RM matrix, for the positive canonic RM expansion and the 
RM matrices for the non-positive canonic RM expansions. The 
generalisation is based on polarity function of the variables 
in the Boolean function. 
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The eigen-structure of the RM matrices would be 
different for different structure of the RM matrices. It 
has been possible to define the number of eigen-values of 
the RM matrices as a function of the total number of variables 
of the Boolean function and the number of variables of these 
that appear complemented. Hence, eventhough the eigen-values 
of two different polarity functions may be same, as the 
number of variables and the number of variables being comple- 
mented being the same for two different polarity functions, 
the eigen— vectors of the RM matrices with these two polarity 
function would not be the same. The RM matrices are 
obtained, as direct product of the core matrices, in sequence 
as per the pattern of the polarity function and hence the 
RM matrices for different polarity functions would be 
different and so would the eigen-vectors be. 

A Boolean function is expressed in binary and if it is 
not possible to obtain the eigen-values in the binary field 
GF(2), then the roots are found in its extension field. The 
characteristic equations of the RM matrices, for the non- 
positive canonic RM expansion, being a product of an irre- 
ducible polynomial of degree 2, it is only natural to expect 
the roots of the characteristic equations to be repetitive, 
being confined to the elements of the field. The method of 
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obtaining the eigen-vectors when repetitive roots occur is 
based on the generalised inverse of the characteristic matrices. 
A specific technique has oeen developed in this thesis to find 
the modal matrix of the RM matrix for the positive-canonic 
RM expansion. 

Eigen-structure of any matrix characterises the 
transform. The eigen-vectors of the RM matrix for the positive 
canonic RM expansion, partitions the vector space C, the 
vector space of all possible sets of the minterm coefficients 
c^s. Vtfith the set of eigen-vector as a subgroup, cosets 
have been formed of this vector space C. By identifying the 
co set-leader and associating the coset-leader with an 
eigen-vector, through the characteristic matrix, the minterm 
coefficient is related to its corresponding positive canonic 
RM expansion coefficients. Similarly, the eigen-values of 
the RM matrix, for the non-positive canonic r*M expansion, 
have been associated with the characteristic matrix and the 
minterm coefficients, to generate the non-positive canonic RM 
expansion coefficients. 

1.3 OUTLINE OF THE.. THESIS: 

In this chapter the aim of this thesis and the common 
terms used in this thesis were explained. The line of 
approach towards the generation of the RM expansion coefficients 
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was also given. The study of the eigen-structure itself was 
considered important because some input vectors remain 
unaltered by the transformation, except by a scaling factor, 
called the eigen-values. These input vectors are called as 
the eigen-vectors and they characterise the transformation 
and these could be put to appropriate use. 

Chapter 2 is a mathematical recaptulation of the pre- 
liminaries necessary in the development of this thesis and 
highlights some techniques for obtaining the eigen-vectors. 

In the next chapter, Chapter 3, the structure of the RM 
matrices have been studied and has been represented as direct 
product of core matrices. Properties of direct product of 
matrices have been applied to find the eigen-values of the 
RM matrices. The direct-product method of finding the eigen- 
values of the RM matrices, for the non-positive canonic RM 
expansion, from the eigen-values of the core matrices have 
been simplified into formulae. Characteristic equations of 
the RM matrices, from the eigen-values so obtained, has also 
been outlined. A simple method of finding the modal matrix 
of the RM matrix, for the positive-canonic RM expansion, has 
also been developed in this chapter. To obtain the modal 
matrix of the RM matrices, for the non-positive canonic RM 
expansion, the generalised inverse of the characteristic. 
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matrix, especially when the roots of the characteristic 
equation are repetitive, has been suitably adopted* The 
RM matrix has also been factorised using Good’s technique* 

In Chapter 4, the minimisation of the Ex-OR design 
has been discussed. Algorithm for generation of the 
positive canonic RM expansion coefficients and non-positive 
canonic RM expansion coefficients, directly from the minterm 
coefficients have been given. These algorithm are based on 
the eigen-vectors of the RM matrices and the eigen-values 
of the RM matrices respectively. Based on Good’s technique 
of matrix factorisation sequential (binary/decimal) genera- 
tion of all RM expansion coefficients has also been done. 

The last chapter summarises the work in this thesis. 
Suggestions for further work in this field has also been 
given in this chapter. 



CHAPTER 2 


MATHEMATICAL PRELIMINARIES 

This chapter briefly outlines the mathematical pre- 
requisites. These include relevant results from the theory 
of groups, rings, fields and matrices. Matrices not only 
simplify the analysis but help in organising computer methods 
for the techniques developed. Here the algebra of matrices 
essential for the thesis has also been included. Standard 
results have been adopted from a number of text books, a 
few of these are [ 9 ], [12], [13] , [14] , [15] , [16], [l7] , [l8], 

[19], [20]. 

2.1 GROUPS, RINGS AND FIELDS: 

2.1.1 Groups: 

A group is a set of elements with one operation and ; 

t 

its inverse, such as addition and its inverse viz. subtraction j 
and multiplication and its inverse viz. division. A group G j 
thus is a set of objects for which an operation is defined 
and for which the following rules hold. Let a,b,c be 
elements of the group. The operation is a single valued 

i 

function of two variables, sayf(a,b) = c, which is | 

denoted as a+ b = c 01 ab = c and called as addition or 
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"multiplication even though it may not oe addition or 
multiplication of the arithmetic of ordinary numbers. 

Rule l (Closure); The operation can be applied to any two 
group elements to give a third group element. 

Rule 2 (Associative Law): For any three elements a,b and c 
of the group, (a+b)+c = a-t(b+c) if the operation is written 
as addition , or a (bc)=(ab)c if the operation is written 
as multiplication. 

Rule 3 : There is an identity element. If the operation 
is addition, the identity elemnt is called zero, '0* and 
if the operation is multiplication, the identity element 
is one, • 1* . 

Rule 4 : Every element of the group has an inverse element. 

If the operation is addition, the inverse element corres- 
ponding to ’a* is '-a' and is defined by the equation as 
a+(-a) = (-a)+a=0. If the operation is multiplication, the 
inverse of 'a* is • a”^ 1 and is defined as aa"”^ = a~^a — 1. 

In addition to the above laws, a group may satisfy the 
commutative law; i.e. a+b = b+a, or if the operation is 
multiplication, ab = ba. Such a group is called an Abelian 
or a Commutative group. 

2 . 1 .2 Rings: 

A ring R is a set of elements for which two operations 
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are defined. One is addition and the other multiplication, 

even though these may not be ordinary addition or multipli- 

f or 

cation of numbers* in order £ R to be a ring, the following 
rules must apply: 

Rule 1 : The set R is an abelian group under addition. 

Rule 2 : (closure). For any two elements a and b of R, the 
product ab is defined and is an element of R. 

Rule 3 (Associative Law): For any three elements a,b, and c 
of R, a(bc)=(ab)c. 

Rule 4 (Distributive Law): For any three elements a,b and c 
of R, a(b+c)=ab+ac and (b+c)a=ba+ca. 

A ring is called commutative if its multiplication operation 
is commutative, i.e. for any two elements a and b ab=ba. 

2.1.3 Fields: A field is a commutative ring with an unit 
element (multiplicative identity) in which every non-zero 
element has a multiplicative inverse. 

2.3 SUBGROUPS AND COSETS: 

A subset of elements of a group G is called a subgroup 
H, if it satisfies all the rules for a group itself. H is a 
subgroup if it satisfies the rules for closure and the rule 
of inverse also applies*. For a set which fulfils the two 
rules, the identity must also be present. 


r!* 
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Suppose the elements of a group G are 9i»9 2 >93»»»» 
and the elements of a subgroup H are h^t^h^..., and consider 
the array formed as follows: 

The first row is the subgroup, with the identity at 
the left and each other element appearing once and only once. 
The first element in the second row is any element not 
appearing in the first row, and the rest of the elements are 
obtained by operating on each subgroup element by this first 
element on the left. Similarly a third, fourth and fifth 


rows are form 

/ 

ad each 

with a 

previously unused group element 

in the first 

column, 

until 

all the group 

elements appear 

somewhere m 

h i 

the array. 

= 1, h 2 

» ^3» • • • 

- h n 

9l \ 

= g r 

g l h 2 

* 91^3* « » * 

' 9 l h n 

«2 h l 

• 

4 | 

= g 2 » 

g 2 h 2 

f g 2 ^3» * • • 

’ g 2 h n 

• 

w m 1 

= V 

g m h 2 

* g m h 3' •** 

> 9 m h n 


The set of elements in a row of this array is called 

a left coset and the element appearing in the first column 

/ 

i'S called the coset leader. Right cosets could be similarly 
formed. Also every element of the group G is in one and 
only one coset of a subgroup H. 
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2.3 GALOIS FIELDS: 

Before defining the Galois fields, associated termi- 
nologies are explained. 

2.3.1 Ideals, Residue Class and Residue Class Rings: 

An ideal I is a subset of elements of a ring R with 
the properties that I is a subgroup of the additive group of 
R and for any element 'a 1 of I and any element ' r* of R> ar 
and ra are in I. Since an ideal is a subgroup, cosets can be 
formed. In this case the cosets are called residue classes 
and all properties of cosets apply to residue classes also. 
The residue classes of a ring with respect to an ideal forms 
a ring and this ring is called the residue class ring. 

2.3.2 Ground Fields, Extension Fields: 

The residue class ring modulo m is a field if and only 
if m is a prime number. These fields are called prime fields, 
or Galois fields of p elements, GF(p). 

If r(x), s(x) and t(x) are polynomials and r(x) s(x) = 
t(x), then it is said that t(x) is divisible by r(x) or that 
r(x) divides t(x), and that r(x) is a factor of t(x). A 
polynomial p(x) of degree n which is not divisible by any 
polynomial of degrees less than n but greater than 0 is 
called irreducible. 
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Let p(x) be the polynomial with coefficients in a 
field F. If p(x) is irreducible in F, that is, if p(x) has 
no factors with coefficients m F, then the algebra of 
polynomials over F modulo p(x) is a field. The field formed 
by taking polynomials over a field F modulo an irreducible 
polynomial p(x) of degree k is called an extension field of 
degree k over F. The original field F is called the ground 
field. 

2.4 GENERALISED INVERSES OF MATRICES: 

The concept of division does not exist in matrix 
algebra, however, an equivalent operation results when the 
matrix is multiplied by an unique matrix, called an 'inverse* 
to obtain the identity matrix, the 'one' of matrix algebra. 

The idea of inverse is that if Ax = b then x = A^b and this 

—1 

A is such that A A = I. This A exists only if matrix 
A is non-singular and is a square matrix. 

This unique inverse is at times denoted as M also, 
as this matrix M satisfies the four conditions given by 
Moore and Penrose and hence at times it is also known as the 
Moore-Penrose inverse. The four conditions are, 

i) AMA = A 

ii) MAM = M 

iii) AM is symmetric and 

iv) MA is symmetric. 
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The matrix M defined by the four Moore-Penrose 
conditions is unique for a given A. But there are many 
matrices G which satisfy just the first Moore-Penrose 
condition* 

AGA — A 

Any matrix G that satisfies this condition is called a 
generalised inverse of A, for A of dimension pxq and G of 
dimension qxp. Thus many matrices G may be found that 
satisfy this condition, the one exception is when A is non- 
singular in which case there is only one G and it is the 
regular inverse G = A~ X =M. The matrix G has numerous other 
names, depending on the various number of conditions it 
satisfies. The M inverse of A is denoted by A"*^ and the 
generalised inverse by A”. 

2.4.1 An Algorithm to Find the Generalised Inverse G: 

Matrix A has dimension pxq and is the leading 
sub-matrix and is non-singular of rank r^. Then a generalised 
inverse may be obtained as 



r, u 

A 12 

f 

r_i 

Vi 

0 1 

A p x q ~ 

i 

! 

i 

A 21 

A 22_ 

is G„ = 

qxp 

u° 

of 


where the null matrices in G have appropriate order to make 
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G q x p, the verification of AGA = a involves using 
A 22 = A 21 A 11 A 12* 

However if is not a non-singular sub-matrix of 

order x, then by elementary operations on the matrix a 
non-singular submatrix of order r may be brought to the 
leading position and then by partitioning the matrix B, 
and applying the same elementary operation on the generalised 
inverse so obtained, the generalised inverse G of matrix A 
may be obtained. Let R and S be the permutation matrices 
such that RAS brings a non-singular matrix of rank r^ to 
the leading position. Thus 



then F, the generalised inverse of B is F = 

s 

and therefore G = SFR, which is a generalised inverse of A. 

2.5 DIRECT PRODUCT OF MATRICES: 

The direct product of two matrices A p x ^ and B m x n 
may be defined as 

A pxq® B raxn 
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and is called as Kronecker product or Zehfuss product. 

The product matrix is partitioned into as many submatrices 
as there are elements of A, each submatrix being B multi- 
plied by an element of A. Therefore the elements of the 
direct product consist of all possible products of •* elements 
of A multiplied by an element of B. It has, in general, 
order pm x qn. 

Direct product, therefore, is an unique manner of 
factorising a matrix defining a transformation and hence 
implementation of the transformation can be accomplished 
with reduced number of arithmetic operations. 

Some useful results of direct product of matrices are 


a) R <ij> “ R <i 0 i 1 l 2 ><3 0 j 1 j 2 > ” 3l 0 3 0 b iPi C i 2 J 2 *** 

where R ~ A0B (&- C ... and are mixed radices. 

b) If a^ and are the eigen-values of A a matrix and 

b^ and b 21 and so on of the others then the eigen-values 
of the resultant matrix a(x)b(x)... is the product of their 
individual eigen-values. 

i.e. a ll b 21* a 21 b ll and a 21 b 21 and S0 on * 


c) If a u , b n 
a 21 b 21 


are simple roots of the characteristic 


equation of the matrices A and B, then the eigen-vectors of 
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of C = A <£ B are a ii b ii and So on. 

a ll b 21 

a 21 b ll 

a 21 b 2l 

2.6 EIGENVALUES AND EIGENVECTORS OF MATRICES: 

For u a vector and y a scalar, the equation 
Au = yu 

may be valid for a given matrix A. The conditions for which 
this equation is valid, other than for u = 0, is that A-yl 
is singular i.e. j A—yl | = 0. The above condition may be 
rewritten as 

(A - yl) u — 0 

A non-null solution for u may be obtained as 
u = [(A - yl)~ (A-yl)-l]z 

using a generalised inverse of (A-yl) and where z is arbitrary. 

The condition |A-yl| = 0 is called the characteristic 
equation of A. For A of order n, the characteristic equation 
is a polynomial equation in y of order n, with n roots 

t 

denoted by y^, V 2 » • * *Y n » some of which may be the same and : 

some of which may be zero. These roots are called latent-roots, | 
characteristic-roots, proper roots, eigenvalues or y-roots. ; 

Corresponding to each y. is a vector u. which satisfies 

Ju J* 
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A ik = y± u i ^ or i = tr*** n a nd these vectors u j_» u 2»*'‘ u n are 
correspondingly called latent-vectors, characteristic-vector 
proper-vector or eigenvector. And the matrix [A-yl] is called 
the characteristic matrix. 

The set of all vectors u, which form the columns of a 
matrix is called as a model matrix. Model matrix may also 
be found using Jord<vn canonical matrix [l2], 

2.6.1 Calculating Eigenvectors: 

Supposing y k i s an eigenvalue of A. It is a solution 
of the characteristic equation |A-ylj = 0. Determining an 
eigenvector corresponding to y k implies finding a non-null u 
to satisfy Au = Y^u, which is equivalent to solving 

(A - y k l)u = 0 

The number of linearly independent solution to the above 
equation are n-r(A-y k I) i.e. the order of the matrix minus 
the rank of the characteristic matrix (corresponding to the 
characteristic value) and these solutions are 

u k = [(A-y k l)~ (A-y k I)-l]z and 

since y k is such that |A-y k l| = 0$ A-y k I is singular and so 
at least one non-null solution always exists. For each y k , 
the characteristic matrix [A-y k l] is processed by elementary 
operations and a matrix B k is obtained to enable partition 
the matrix as, 



7 



such that the submatrix is non-singular with the same 
rank as the characteristic matrix [A-y^l] and then obtain 
the generalised inverse of the characteristic matrix as 



and so obtain the solution 



A number of eigen-vectors are found from the 
general solution for each eigenvalue y k of the matrix 
including if it occurs. Since y k is a solution to 

the characteristic equation it may be a solution more than 
once, in which case it is called a multiple root. If y^ is 
a solution only once it is called as a simple root. When 
y k is a simple root the rank of the characteristic matrix 
r(A-Ykl) 38 n-1. Therefore there is just only one LIN eigen- 
vector corresponding to y^. This is reflected in w as w 
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becomes a scalar. When y^ is a solution more than once 
it is called as multiple roots and the number of times that 
it is a solution is called its multiplicity. For each 
multiple eigenvalue y^, with multiplicity m^, the rank of 
A-y^I must be ascertained as this is the value that plays 
an important role in the determination of eigen-vectors. 


2.7 NULL-SPACE AND DIMENSION OF NULL-SPaCE OF MATRICES: 

In Sec. 1.1, the linear transformation was identified 
with an n x n matrix R. T was the linear transformation of 
a vector space C into a vector space a. The null-space of a 
linear transformation T is the set of all vectors such that 
= o, and the null-space of a matrix R is the set of all 
row matrices X which satisfy the homogenous linear equation 
XR = 0. Alternatively, the null-space of the row space of a 
matrix is called the null-space of the matrix. A vector is 
in the null-space of a matrix if it is orthogonal to each row 
of the matrix. The dimension of the null-space of a given 
matrix is called as nullity of the matrix. The rank and 
nullity of an n x n matrix is related by thv following 
equation. 

Order of the matrix - rank of th ©matrix = nullity. 



CHAPTER 3 


EIGEN-STRUCTURE OF REED-MULLER MATRICES 

Methods to obtain the positive canonic RM expansion 
coefficients from the mmterm coefficients (obtained from the 
truth-table entries of a Boolean function) have been 
developed [7] , [lO] , [21] , [22] , [23] . In the methods developed 
in [7],[l0] other non-positive canonic RM expansion coefficients 
have been obtained from the positive canonic RM expansion 
coefficients. In the geometric mapping method [7], the map 
folding operations become difficult when the number of 
variables in the Boolean function increases. Moreover, if 
this difficulty is surmounted then in any map method there 
is no freedom to represent the variables than the available 
3 dimensions, with utmost four variables represented in each 
dimension. The limitations of the b^ maps have been overcome 
m the computer method [lO]. 

In this chapter the RM matrix structures, for different 
polarity functions have been obtained as a direct-product of 
two kinds of core matrices. Techniques have also been 
developed, or suitably adopted, to determine the eigen-structure, 

3.1 RM MATRICES: 

-> ^ 1 

In Chapter 1 the method of obtaining any non-positive > 
canonic RM expansion from any other expansions by algebraic 
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substitutions and simplification [24] was illustrated by 
(1*8), (1.9), (1.10) and (l.ll). For this example the 
transformation, in the matrix form wo uld be 

01000000 
11000000 
01010000 
11110000 
R= 01000100 
11001100 
01010101 
11111111 


This is the transformation from tht, mxnterm coefficients, 
obtained from the truth-table entries of the Boolean function, 
to the non-positive canonic EM expansion coefficients. This 
particular transformation pertains to the polarity function 
110 of the 3 variable Boolean function. The RM matrix, for 

I 

the positive canonic EM expansion, was also obtained iteratively) 

I 

(1.14) which may also be obtained as direct product of, what j 

i 

has been termed as the core matrix, for m' = 0 in (1.14), i.e. j 
the 2x2 matrix [f . To obtain the RM matrix, for the , 

non-positive canonic RM expansions, suitable core matrix, { 

( 

to be called as core matrix for a complemented variable, j 

would have to be appropriately multiplied with the other core j 


I 
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matrices; by multiplication is meant the direct product 
as per some rules. 

The eigen-values and eigen-vectors of the Reed-Muller 
matrices help in obtaining the positive canonic and non- 
positive canonic expansion coefficients and therefore, the 
eigen-structure of the RM matrices have been developed in 
this chapter subsequently. The results obtained in this 
chapter have been used m the next chapter, m which tne 
algorithms to obtain the expansion coefficients have been 
developed. As the Reed-Muller matrices have been developed 
as direct product of the core matrices, the eigen-values 
too have been obtained as direct product of the eigen-values 
of the core matrices. As brought out earlier the eigen-vectors 
for repetitive roots (eigen-values) of the characteristic 
matrix of RM matrices, have to be uniquely obtained by 

i 

employing the generalised inverses of the characteristic ; 

matrix or using the concept of null-space of the characteristic i 
matrix. I 

3.2 THE CORE MATRICES: 

I 

In the previous section the core matrix for the > 

1 positive* variable was identified. Like the RM matrix, 

for the positive canonic RM expansion, this core matrix is also j 

involutary, i.e. the matrix is its own inverse. Similarly t 
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to obtain the RM matrices, for non-positive canonic RM 
expansions, a core matrix to represent the complemented 
variable is identified as that matrix which results when 
the core matrix for a positive variable is multipled by the 
counter identity matrix J. This is represented as 


= b x .j = “j [“ J] 


[° i ] 

L 1 1 J 


and may be called as the core matrix for a complemented 
variable* 

The characteristic equations of the two core matrices 
are respectively, 


(y + 1)^ = 0 (3.2*2) 

( y 2 +Y+l) = 0 (3*2.3) 


The roots of the characteristic equation of the core matrix 
is y=l and y=l. The repetition of a root is also termed as 
its multiplicity and is denoted by m k . Therefore the 
characteristic roots of the core matrix may be written as 
y=l, m^=2. The characteristic equation of the core matrix 
for a complemented variable is an irreducible polynomial of 
degree 2, Therefore, the roots are not in the ground field, 
GF(2) and are to be found in the first extension field viz. 
(GF2^) and the roots would be from the elements of the 
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extension field, viz. 0,1, a and 1+a, with a^+a-KL = 0 being 
the irreducible polynomial. Hence, the ch: racteristic roots 
of the core matrix of a complemented variable are y=a and 
y-l+a, with multiplicity one each. It is also observed that 
the irreducible polynomial is of degree 2 and therefore the 
characteristic roots always occur as a set. 

3.3 OBTAINING RM MATRICES AS DIRECT -PRODUCT OF CORE MATRICES: 


A given m variable Boolean function may be associated 
to the corresponding RM expansion coefficients through trans- 
forms, called as the RM matrices. It was also described in 
Chapter 1, that the RM matrices may be classified as the RM 
matrix for the positive canonic RM expansion or the CFRM 
expansion and as the RM matrices that represent the other 
non— positive canonic RM expansions, called as the generalised 
Reed-Muller expansions. In this section, the method for 
obtaining the two kinds of Reed-Muller matrices, as direct- 
product of the two kinds of core matrices, is developed. 

3.3.1 RM Matrix for Positive Canonic RM Expansion: 

The RM matrix for the positive-canonic RM expansion, 
in which all the m variables of the Boolean function are 
in the positive form, is obtained as direct product of the 
. The m variable 3oolean function is said to 


core matrices 
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have all positive polarity function and hence each of the 
m variable of the Boolean function is represented by a core 
matrix* Therefore, the Reed-Muller matrix, so obtained 

is of order n, n=2 . For the Boolean function of m variables, 
this may be written as 

R - f-i ?3 (?) <?> .... rx) [* °] m times (3.3.1) 

The RM matrix so obtained is also mvolutary. 

3.3.2 RM Matrix for Non-positive Canonic RM Expansions: 

In the non-positive canonic expansion, one or more 
variables of the Boolean function apoear in the complemented 
form, but the same variable does not appear m both the forms 
in the expansion. The RM matrix corresponding to these 
expansions would, therefore, be dependent on the polarity 
function of the Boolean function and is obtained as direct- 
product of both kinds of core matrices. Each variable in the 
function which has positive polarity is represented by the 
core matrix and the variable that is complemented by the core 
matrix of a complemented variable. For the Boolean function 
of m variables, which has the polarity function 1101001... 1, 
the RM matrix is obtained as direct product of the appropriate 
core matrices in that order. 
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R = [£?]© [ 


1 0 
1 1 


j«i? i 


] ® t 


1 o 

l l 


j ® [j 


l 

l 


^ [° l] ® [ 1 X 


] 


(x) ... 0 [| ^] it times (3.3.2) 


Therefore, if there are t variables complemented m the m 
variable functions, there are (m-t) core matrices and t 
core matrices, for a complemented variable. The RM matrices 
so obtained is unlike the RM matrix for the positive canonic 
RM expansions, i.e, is not involutary. Here too the order of 
the RM matrix is n, n=2 m . 

3.4 GOOD ' S FACTORISATION OF RM MATRICES; 

The RM matrices have a high degree of redundency and 
if the redundency in the definition of the matrix trans- 
formation is eliminated by matrix factorisation then a more 
efficient means of implementation is possible. Such a 
technique was employed in the development of the fast Fourier 
transform (FFT) and fast Hadamaxd transform (FHT) and to a 
large class of fast transformations [14], 

A class of matrices formed by Kronecker product 

operation may be considered as the square submatrices of 

order p, with entries in , where i and j range from 0 to 

I > 1 ( J 

p-1, and the first index r representing the class of entries 
corresponding to a particular dimension in the Kronecker 
product operation. In general 
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m 


m„ ~ i 

r,o,o r,o,l 
r,l,o m r,l,l 


M r = 


and 


m 


m_ , „ m_ , , . . . m, 


r,o,p-l 

x r l,p-l 


m r p-l o m r, p-1,1* * * m r, p-1, p-1 


H 1 = M o 


H = M 1 ©H 1 


( 3 . 4 . 1 ) 


( 3 . 4 . 2 ) 


Hq = M q-1® H q-1 


Thus 


H. 


m. 


H, 


l o ^q-1 m q-l,o,p-l m“1 


H„ 


m q-l}l,o H q-1 *’* r V-l>l»P“ 1 “ q_1 


m 


q— 1 , p- 1 , o H q-1 


m. 


H. 


q-1, p-1, p-1 


( 3 . 4 . 3 ) 


where 


H is a p q ^P q matrix 


When operating with Kronecker matrices within a 
computer, it becomes 'desrrable to store a representation 
(algorithm) of the entries of the product matrix rather than 
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the matrix itself. The location m the matrix may be described 
by their lexicographic sequence representation# A given 
index of matrix may thus be represented by q digits each 
of which may take the value 0 to p-1. Representing the 
horizontal index by v and the vertical index by x, the names 
of the rows and columns in lexicographic sequence for 
matrix for p=3 is 


v 


00 01 02 10 11 12 20 21 22 


00 

X 

01 
02 
10 


HjT 11 
12 
20 
21 
22 




H 2 (*,v) 


(3.4.4) 


Representing the v and x variables in the lexicographic number 
system, modulo p Requires q digits to allow v and x to range 


over 0 to p * . 

Therefore v and x 

may be described 

by 

v = v n-l 

V2 • ’ • v X v o V i 

Q,l,***> P~1 

(3 

x * x n-l 

x n-2 •" x l x o x i 

0,1, • • • , p— 1 



f 
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Using such notation allows the entries of the p x p matrix 
(3.4,2) to be described by the equation 


H 1 (x,v) 


P-1 

7t 

i=o 


P-1 

it 


(x 0 -i)d(v o -j) 


(3.4.6) 


where d(a-b) is the delta function which takes on the. value 
1 whenever a=b and 0 otherwise. The representation of 
(3.4.6) may be interpreted as multiplying oil entries of the 
core matrix, M Q (3.4.2), together and noting that all but 
one entry would be raised to the Oth power. The- entries of 
the p by p matrix, lUt (3.4.4) may now be represented as 


^2 ( ^ = 


P-1 P-1 

71 71 

i=o j-o 


m 


d(x 1 -i)(v 
l»i» j 


1 


j) 


P-1 

% 

i=o 


P-1 

71 

J=0 


d(x -i)d(v -j) 

w U O 

o, i , 3 

(3.4.7) 


where the exponents determine the correct product of entries 
for a given v and x. In general, the entries for H q may be 
represented as. 


H q ( x , y) 


q-1 P-1 P-1 d(x -i)d(v r -j) 
n 7i % m . 
r=o i=o j=o I,x ’ 3 


(3.4.8) 


following tht recursive notation of (3.4.6) and (3.4.7). 
Representation of the rows or columns of a Kronecker matrix 
in the form of (3.4.8) allows the generation of any single 
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element, column or row of the matrix without storage of the 
entire matrix array. The operation indicated in (3.4.8) 
may be visualised as multiplying all entries of the sub- 
matrices, M r , forming the Kronecker Product, together and 
letting the exponent operation allow only the proper entries to 
be raised to power 1 while all others are raised to the power 

0. For the Kronecker product of identical matrices 

n-1 

p-l p-l l d(x r -i)d(v r -j) 

H ( x, v) - it n m* 9 (3.4.9) 

4 i=o j=o * J 

A major observation was made about the matrix facto- 
risation i.e. if highly redundant matrices could be factored 
into a product of matrices with few non-zero entries, then a 
fewer number of operations would be necessary for transfor- 
mation implementation. A technique, attributable to Good 
[26], of matrix factorisation^may be used to decompose the 

class of Kronecker matrices described by (3.4.8) or (3.4.9). 
the 

Thus for Hq Kronecker matrix of (3.4.8), there exists q 

matrices, each of dimension p^, such that when multiplied 

together will equal H . These matrices may be described as 

9 
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In this matrix there are p 0 **^ non-zero entries and only 
non-xedundant elements. Then 


H n “ ^- G n-1^ ^ G n-2^ * * * f- G J[ G 0 J 


(3.4.11) 


and for the Krocker 


H n = 


product of identical matrices 
[G] n (3.4.12) 


The RM matrices obtained earlier (3.3.1) and (3.3.2) 
may be factorised respectively as, 
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whore. 


R 


10 


10 


10 


10 


11 


11 


11 


11 


•11 


n columns 


n/2 rows 


n/2 rows 


and, 


R = 


( 3 . 4 . 14 s ; 


— _ 

2 

— 

- 1 

10 


01 


10 

11 


11 


11 

10 


01 


10 

11 


11 


11 

10 


01 


10 

11 


11 

* 


11 

* 

* 

* 

• 

I 

I • 

9 


• 

• 

10 

1 

1 

01 


10 

11 

| 

1! 


11 

— * 

L 

— 




. as per 
the 

polarity 

function 


( 3 . 4 . 15 ) 


= R^ R^ . R~ R^ Rj *** m ^^ mes 
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where R. =* R and R. = R- as per the polarity function, and 

X a X X 

01 

rji 

01 

01 

R — — 

x 11 

11 

11 

♦ 

11 

both R and R- obtained by rearranging the rows - e -f (3.4.10). 

X X 

3.5 EIGENVALUES OF THE REED-MULLER MATRICES; 

The eigenvalues or the characteristic roots of a 
matrix are the solutions to the equation, 

Ru = yu (3.5.1) 

for R a matrix and u, a vector and y, a scalar. There are 
a number of scalars and vectors that satisfies this equation; 
the scalars are called as eigenvalues and the set of eigen- 
values is are termed as the spectrum of R. The eigenvalues 
are also obtained by solving thecharacteristxc equation 

( R - yl) = 0 (3.5.2) 

In the previous section the RM matrices were obtained as direct 
product of the core matrices and were classified as RM matrix 


n/2 rows 

(3.4.16) 

n/2 rows 
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fox the positive canonic RM expansion and RM matrices fox 
the non— po si tive canonic RM expansion. The polarity 
function of the variables of the Boolean function determines 
the kinds of roots and their multiplicity. In this section 
the eigenvalues are not found as solutions to the chara- 
cteristic matrices of the RM matrices (3.5.2) but the property 
of direct-product is exploited. 

3.5.1 Eigen-value of the RM Matrix, for Positive Canonic 
RM expansion: 

The RM matrix for the CFRM expansion was obtained as 
direct product of the core matrix [^] (3.3.1); the order of 
the matrix being n=2 m . Therefore there are utmost 2 m 
eigenvalues. The eigenvalues of the core matrix are y=l, 
with multiplicity m^=2 and therefore, the eigenvalues of the 
RM matrix for the CFRM expansion are y=l, m^n. And, 
therefore may be written as 

Spectrum of R = [^] [^] ® *•* [j] » y=l,m 1 =2 m (3.5.3) 

3.5.2 Eigen-values of RM Matrices, for Non-positive 
Canonic RM expansions: 

The eigenvalues of the RM matrices, for the non-positive 
canonic RM expansion may be obtained as the direct product of 
the eigenvalues of the core matrices. The characteristic 
equation of the core matrix for a complemented variables being 
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an irreducible polynomial of degree 2, the roots are found 
in the extension field only and the roots occur as a set 
viz. y — a and • =l+a. The number of eigenvalues of the RM 
matrices is dependent on the polarity function i.e. on the 
number of variables of the function and the number of 
variables complemented. Hence, the number of eigen-values 
of two RM matrices may be the same but the pattern of the 
polarity function being different, their eigenvectors would 
be different. 

For the polarity functions in which there is only one 
variable complemented, all roots of the characteristic 
equations for the RM matrices are to be found in the extension 
field alone. This is unique only to the case mentioned, 
because irrespective of the number of variables, if there is 
only one variable complemented of these, then the eigen-values 
of the core matrix of a complemented variable manifests 
the entire operations of the direct product and therefore 
the roots are found only in the extension field. 

In general it may be written as 

R = [J°] ® ... ® *>• 

Spectrum of R - 

[Jj <$) ... [J] & Ci«] ® tj] "i times 
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l(2 m “ 2 )' 

l(2 m - 2 ) 


to 


a 

1+a 


a {2*" 1 ) 


1-Hx(2 m *’ 1 ) 


i.e. Y 
Y 


a, m ff = 2 


,m-l 


1-kx , m 1+a * 2 


,m-l 


Cases where more than one variable occurs complemented 
the eigenvalues y=l, y=a and y=l+o are obtained and satisfy 
the condition 


<“l + % + “i+a = 2 


,m 


<3 . 5 . 4 ) 


as utmost 2 m eigenvalues may be had for the RM matrix of order 
n , n=2 ra . As a result of the polarity function, the first three 
variables in this category may occur as 

l\i 8> c 1*3 ® ^ - 


1(4 

«(2) v 
l+a(2) 


or 


& [ i ] 


l+a* 


or 


IJU ® & ® tiW 


8»L 

l+a(2) 


1(4 

a(2) 

l+<x(2) 


is being indicated in parenthesis) 


(for ease, the multiplicity 
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Now the variables following may be true or complemented, and 
would occur as per the polarity functions, if the variable is 
true then 


a(2 
l+cc(2) 



1 ( 8 ) 

a(4) 

l+a(4) 


( 3 . 5 . 5 ) 


i.e 


m 1 ' = 2.11^, 


K = 2.ra a and ra' 1+a = 2.ra 1+a 


(where the prime on the multiplicity indicates the multiplicity 
after the direct product operation) and if the variable is 
complemented, then 


1(4 
a( 2 
1+a ( 2 ) 




1(4) 
a{6) 
l+a( 6) 


(3.5.6) 


1 . G . 


ra i * V™!**’ 


m f =m, +m . 
a 1 a* 


m 


1+a 




It was observed that the eigenvalues, though a function 
of the polarity function, depended more on the number of 
variables complemented; the odd number of complemented variables 
demonstrating a common feature amongst them and the even 
number of complemented variables demonstrating another 
common feature amongst them. 

Starting with two complemented variables, the resulting 
eigenvalues would be 



50 


[ 




a 1 
l+a J 


1 ( 2 ) 

a(l 

l+a(l)_ 


(3.5.7) 


then with other complemented variables following this, the 

rule given in (3.4.6L ) could be applied. It was observed 

that when even number of variables are complemented 

m l ~ m a + ^ “ m l+a +1 and when ocld number of variables are 

complemented then = m a -l = m 1+a -l, which results as an 

outcome of the operation over the irreducible polynomial 
2 

y +y-**l=0. The eigenvalues of ten complemented variables 
are considered before generalisation. 


The eigenvalues of t complemented variables, obtained 
as a direct product of the eigenvalues of the core matrices 
for complemented variables, is tabulated to illustrate the 
observation that the eigenvalues are indeed a factor of t being 
odd or even and that the rules 

m l = m a +1 = m l+a +1 ( for 1 even ) 
and m^ = m a -l = (f° r * ocicl ) holds. 

The table shows the multiplicity of each root. 
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Table 3.5.1 


No. of 
varia- 
bles 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

• 

The 

eigen- 

values 

a 

‘ 1+a 

a 

1+a 

a 

1+a 

a 

1+a 

a 

1+0 

a 

1+a 

a 

1+a 

a 

1+a 

a 

1+a 

a 

1+a 


1 

B 

fl 

2 | 

6 

■p 

it 

■■ 

22 

42 

86 

170 

| 342 


a 

B 

s 

3 

5 

u 

21 

43 

85 

171 

: 341 


' 1+a’ 

B 

B 

3 

5 

immmmm mm 

21 

43 

85 

171 

341 



Now generalising (3.5.6), two sets of forrrulae on<- each for 
t odd and t even may be written as 


l ) t odd: 


Y = 


Y = 


. 2 t -2 
1 m x = -3-- 

2 t+1 +2 

a m ~ = ““ 


a 


2 t+1 +2 


y = l+o m 1+a = —5 


(3.5.8) 


ii) t even: 


2 t +2 


Y = 1 m x - ™“3 

Y = a m o; “ 


2 t+1 +2 


Y = 1+a m 


2 t+1 -2 


1+a 


** 




(3.5.9) 


8761)8 


•wai > 
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If t = m, i.e. 
the condition (3.4.4) 

m l +m a +m l-Kx = V 2 + 

2* +2 . 
= ~T“ + 


all variables are complemented then 
must also hold i.e. 

2 t+1 +2 . 2 t+1 +2 rt t „m 
6 + ~ * L ~ z 

2^=2 + 2 !!i =2 „ 2 t . 2 m 

D O 


Now if these t variables were complemented out of 
the total m variables of the function, the remaining (m-t) 
variables being true, the eigenvalues of the RM matrices is 
obtained by multiplying the eigenvalues obtained for t 
variables by a factor of Eqn. (3.5.3) refers. 

Therefore the multiplicity of the eigenvalues of RM matrices 
for an m variable function, with t variables complemented as a 
function of m and t are obtained as 


i) t odd : 


Y 


, oU-t) 2-2 

» 1 mj_ = 2 V “5 — 


(3.5.10) 


0 (m-t ) 2 t+1 ±2 

y = a m a = 2 V —3 


_ _ «(m-t) 2 t - fl ±2 

y * 1+a m 1+a - 2 6 


ii) t even ; 

_ _ ,(m-l) 2^+2 

m l “ 2 3 


Y = 1 


(3.5.11) 
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Y = a m a =2 


(ra-t) 2 t+1 -h2 


Y = 1+a m 1+a = 2 


_ 0 (m-t) 2 t+1 +2 


and these sets of formulae also satisfy the condition (3.5.4) 


m-j+rn + m. . = 2 m , i.e. 

1 a l+a 


t+1. , r» o^+l 


e n(m-t) /2_^2 2 OTX +2 , 2 ^+2 , _ 

«e.,z v o +• r ) - 




and 2 (mrt) ( 2l±2 + 2^2 + 2^2 ) , 2 m . 


Illustrative Example 3.5.1 

The number of variable m = 3 
with polarity x^ *2 X 1 

The corresponding Reed Muller matrix R = C^J ® ^11^ 


00 

10 

00 

00 

00 

11 

00 

00 

10 

10 

00 

00 

11 

11 

00 

00 

00 

10 

00 

10 

00 

11 

CO 

11 

10 

10 

10 

10 

11 

11 

11 

11 
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the eigenvalues are obtained as a direct product of 

j_ (*) i+, a © i anc * the s P®ctrum of R is (a ,a, l+a,l+a,a,a, 
1+a, 1+a ). 


For the same number of variables if the polarity now is 
made as * 3 * 2 x i 


the Reed Muller matrix R = [^] (x) [^] 0 [^] 


00 

00 

10 

00 

00 

00 

11 

00 

00 

00 

10 

10 

00 

00 

11 

11 

10 

00 

10 

00 

11 

00 

11 

00 

10 

10 

10 

10 

11 

11 

11 

11 


The eigenvalues are obtained as j^ a (x) i ©i a direct 
product and the spectrum of R is (a, a, a, a, 1+a, 1+a, 1+a, 1+a) . 


It is thus seen that when there is only one variable 
being complemented, then all roots are m the extension field 
only and irrespective of the pattern of the polarity 
function, the eigenvalues are the same. 

Illustrative Example 3.5.2? 

Consider a 4 variable switching function 
with polarity x 4 x 3 x 2 ^i» 
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the corresponding Reed-Muller matrix R = 
[^]®[^] a direct product 


The eigenvalues are obtained as ( ^ (x) ^ (x) J @ ^ 

the spectrum of R being (a, a, a, a, a, a, a, a, 1+a, 1+a, 1+a, 1+a, 
1+a , 1+a, 1+a , 1+a ) and for the same number of variables, with 
the polarities now assigned as x 4 * 3 x 2 x i* 


The RM matrix becomes R = [^J (x) [^] (T) C^](x) [^] 


The eigenvalues are obtained as a direct product of 

i©A(5»i©i 


The spectrum of R being (a, a, a, a, 1+a, 1+a, 1+a, 1+a, a, a, a, a, 
1+a, 1+a, 1+a, 1+a) . 

It is once again observed that when there is only one 
variable being complemented, then all roots are m the 
extension field only and irrespective of the pattern of the 
polarity function, the eigenvalues are the same. 

Hence for 

m - 3, y=a, m a =4 , and y=l+a, m 1+a =4 and 
m = 4, y=a, m a =8, and y=l+a, m 1+a =8 . 


i 
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Illustrative Example 3.5.3: 

A 7 variable function with five variables complemented 
with the polarity function 1000 100 is considered. 

The function is 

( Xy, x^) and 

the eigenvalues are obtained as direct product of the eigen- 
values of the core matrices as 

c}3 © [£,] © t^]® hM £}]© tiS.] 


the number of 
eigenvalues 
after each 
direct-product 
operation 


1(2) 

a(2) 

1(4) 

1(4) 

1(8) 

1(24) 

1(40) 


o 2 (2) 

«t(fc) 

a(6) 

a(12) 

a(20) 

a(44) 




a 2 (6‘ ) j 

a 2 ( 12 ) ] 

| i 2 (2G) | 

j a 2 ( 44 ) 

1 _ 1 


and by formulae, t being odd 

for y = 1, mj, =* 2^~^ 3"”^ = ~ 40 

for y=a , m a = 2 (7-5) = 4.11 = 44 

for y . 1«, m 1+a - 2< 7 ' 5) ^ - «• 

Therefore m 1 =40, m a =44 and m 1+a = 44 and E m ± = 2 m . 

The same were obtained as direct product of eigenvalues 
of the core matrices also. 



57 


If for the same function the polarity is 1010100 
then m=7, t=4. 

for y=l = 2 ^ 7 " 4 ^ . £-+2 _ 3.5 48 

for y~ a ™ a = 2^ 7 ” 4 ^ . = 8.5 = 40 

for Y~l +a m 1+a =2^ 7 " 4 ^ . = 8.5 = 40 

and Z - 2 m 

3.6 CHARACTERISTIC EQUATION: 

The condition for a non-null solution to the equation 

(3.5.1) 

Ru = yu 

for R, a matrix and u, a vector and y, a scalar is that 
R - yl is singular i.e. |R - yl| = 0 . This condition is 
called as the characteristic equation of the matrix R. For 
the matrix R of order 2 n, the equation is a polynomial 
equation in y of degree n. The roots of this characteristic 
equations are called as the characteristic roots or the 
eigenvalue ( Secs. 2.6. and 3.5 refer). 

For the RM matrices, which are obtained as a direct- 
product of core matrices, the eigenvalues have also been 
obtained as direct product of the eigenvalues of the core 
matrices* These^very eigenvalues being defined by the 
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characteristic equation, it is feasible to’ construct the 

characteristic equation fwmthe eigenvalues, if Y 2 * * *Y k 

are the different eigenvalues with multiplicity in. ,in . . .m , 

Y 1 y 2 Y k 

then the characteristic equation may be written as 
m Y m Y m Y 

(Y“Y X ) 1 (y-Y 2 ) 2 ••• (Y-Y k ) k = 0 (3. ‘6,1) 

And in the case of RM matrices, whose roots are over the 

extension field GF(2^), with the irreducible polynomial 

y^+'Y+l * 0, and all operation being modulo 2, the equation 

may be written as 

m. o m i o 

(y-fl) 1 (y 2 +Y+1) * = (Y+l) 1 (Y +Y+1) = 0 

(3.6.2) 

and these equations may'be specifically written for the 

m variables, of which t are complemented as 

m, m 

(Y+l) 1 = 0 for t = 0 = 2 

( y 2 +Y+l) a -0 for t = 6 m a = 2 m “’ 1 
m, o m i ± l 

and (y+l) *(Y 2 +Y+1) 1 = 0 for 1 < 't 1 m 

and +/- accordingly as t is odd/even respectively. 

3.7 EIGEN -VECTORS OF THE REED-MULLER MATRICES: 

The roots of the characteristic equations of the 
Reed-Muller matrices, being limited to the three kinds of 
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roots that may be had and hence their occurance with 
multiplicity greater than 1, suggests that it may not be 
feasible to uniquely associate an eigenvector to each root. 
Eigen-vectors corresponding to multiple roots may be found 
using the Jordan canonical form of matrices [13]. A more 
specific method of calculating the eigen-vectors, when the 
characteristic equation has multiple roots is by the use of 
the generalised inverse of the characteristic inmatrix of 
the RM matrix R. 

The condition for which non-null solutions of the 
vector u, may be obtained for Ru = Yu is that (R-yl) is 
singular. For singular matrices, only the first of the four 
Moore-Penrose conditions for the unique inverse may be 
satisfied and hence only the generalised inverse may be 
obtained. 

Using this generalised inverses of the characteristic 
matrix and an arbitrary matrix z, the non— null solution of 
u may be obtained as 

u = [(R©Yl)“ (R®yl)@l3 z (3.7.1) 

Further, the number of such solutions that (R-yI)u=0 may 
have are determined by the order and rank of the chara- 
cteristics matrix for that root, which may be written as 

n-r (R®yD (3.7.2) 


L 
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i.e. these may linearly independent eigenvectors may be 
obtained and with this as the basis, the entire span of 
the eigen-vectors may be had. 


The method of obtaining the modal matrix, i.e. the 
set of all non-null u as the columns, the algorithm given in 
Sec* 2,6 may be used. The characteristic matrix (R-yl) may be 
partitioned as _ _ 


R 0Y k I =B k = 




(3.7.3) 


for each eigenvalues y k , where R fc is non-singular with the 
same rank as the characteristic matrix (R-y^I). And, if R k 
does not satisfy this criteria, then by operating with the 
permutation matrices on the characteristic matrix, R k may 
be made to fullfil this condition (Sec. 2.8 refers). Then 
the generalised inverse obtained for the characteristic.' 
matrix is 


(R®Y k l)~ = 


R 


0 


-1 


0 

0 


(3.7.4) 


Further, considering the arbitrary vector Z and partitioning 
it as 

, Z’ = [— V, Jf'W]' (3.7*5) 

and substituting in (3.7.1), a number of non-null solution 
governed by (3.7.2) for arbitrary w of order n-r(R-y k l), 


is obtained as 
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u k = 


+R k C k 


W 


( 3 . 7 . 6 ) 


The modal matrix so obtained, may be modified to fit 
the field over which all operations are being implemented 
to enable obtain the same with ease. The arbitrary matrix 
w in the modal matrix (3.7.6) may be qualified in an unique 
manner to ensure that the columns are indeed unique thus 
giving the n-r(K-y k l) unique vectors, being the set of 
linearly independent eigen-vectors. The property of the 
arbitrary matrix w, therefore woulct have to be full-column 
ranked. Such a matrix in GF(2) would be an identity matrix. 
And, hence the mode-1 matrix for the Reed-Muller matrices 
may be written as 



(F^-htcgd) 


(Rk'b( c k ) 

\ = 

I 

— 

I 



- 

— — ~ 


The order of the identity matrix being n-r(R-yj c l) end where 

R " ’ is the non-singular submatrix of the characteristic 
K 

matrix of rank, rOv-y^l) and the matrix C^, another 
submatrix of the characteristic matrix whose rows are the 
rows of the R fe matrix extended and the columns are the 
columns remaining from the characteristic matrix after 
the matrix R^ has been obtained. 
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This method of obtaining, the modal matrix may be 
applied to both simple and multiple roots and also where 
the roots are m the extension fields. 

Illustrative Example 3.7.1: 


R = 


For a 2 variables Boolean function with polarity function 11 
the FiM matrix for the positive canonic KM exoansion is 

10 0 0 
110 0 
10 10 
[1 1 1 1 

and the characteristic matrix is 

0 0 0 0 
10 0 0 
10 0 0 
1110 


(H-yD 


1 


the order and rank of the characteristic matrix are 4 and 2 
respectively* 


B, 


R. 


D 


'k 



*10 00 ” 


11 10 


00 00 


00 00 
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and R k being involutory I k "’ 1 = R k 

rlOi rOO j _ rOO t 
L 11 J L 10 J “ L 1C J 


and R k ""^C k 


and matrix w, which is arbitrary of order n-r(i - yl) could 
be choosen as the identity matrix of this order, viz. 

- 10 - 


w = [£] 


r 


— 

R k “ lc k w 


rOO-i rlQn 
L 10 J L 01 J 

w 


10 



01 



— ' .1 nr 


and thus 

0 0 
1 0 
1 0 
0 1 

illustrative Example 3.7.2: 

Considering the 2 variable Boolean - function with 
polarity function 10. The RM matrix is 


R = [£] © [?}] - 


0100 

1100 

0101 

1111 


The characteristic matrix 
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(R~Yl) 


y 1 0 0 

1 Y 2 0 0 

0 1 Y 1 

2 

1 1 1 Y 


the rank of the characteristic matrix is 
r(R-Yl) - 3 

Therefore there are n-r(R-Yl) LIN EVS which is 4-3-1. 
Partitioning the characteristic matrix, after permuting 
we obtain 


B = 


1 Y 
0 1 
1 1 
Y 1 


0 0 

Y 1 

r 

1 ^ 

0 0 


for the root y = a 


r 


R, 


a a a 
1 o? 1 


a 2 1 


a 


and R k = 


T 1 

ro 

o 

1 


•m 

0 

0 1 Y 

and = 

1 



2 

1 1 

i 

J 

Y 


with w = [l] 


Y=a 


and hence the model matrix for Y " a 


— * 2 -I 


0 

a a a 0 



1 a 2 1 1 


0 

2 

— 

2 

CN 

i 

CN 


a 

a l a 



U) 

ll 

1 

j 


_ — 


L. 


which indeed is an 
eigenvector. 
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for y = (l+a) the eigen vetted may be found as 
0 
0 
a 
1 

3.7.1 Eigen-vectors of RM Matrix, for the Positive-Canonic 
KM Expansion: 

However, another algorithm for finding the set of LIN 
eigenvector of the RM matrix which represent the positive 
canonic Reed-Muller expansion has been developed. This 
algorithm is based on the concept of null-space of a matrix. 

From the definition of eigenvalues and eigenvectors 
given in Section 2.6, viz. 

Hu - yu for u a vector and y scalar 
it may be written as 

(K - y^) u = 0 

and all vectors u that satisfy this condition forms the 
modal matrices, with the vectors being the columns of this 
matrix. The null space of a matrix may be defined as the 
set of all column matrices x which satisfies the homogenous 
linear equations Dx = 0. An easier approach to obtain the 
null soace of the matrix is by making the matrix B xull 
row— ranked and finding a matrix C which satisfies B C — I 
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and taking those columns that satisfy Bx=0, the null-space 
may be obtained. For a matrix B = (R-yl) this may be written 

as 

(3.7.8) 


B i 

D’ *[ B ^] and C - [C^] 


such that 


B l C i© B 2 C 2 = I <3 ' 7 ‘ 9) 

or equivalently the same equation may be written as, 


B 11 B 12 


C 11 C 12 


”l 0 

B 21 B 2 7- 


C 21 C 2 a 

j 

0 I 


(3.7.10) 


or those C 12 and C 22 that satisfy the conditions 

_ y^S q p =5 0 (3#7*ll) 

B 11* C 12© B 12 C 22 

and B 21* C l2® B 22 C 22 = 1 


q the 

such that c 2 = [ c 12 ] L modal raatrix 

of the matrix rt and being the null-space of the character!- 
stic matrix (R-yl)* 

It is observed that for C 12 = Q and C 22 = I of 
appropriate order, the modal matrix may be formed. The 
choice of W and I also depends on the choree of and B 22 

which satisfy the conditions 


(3.7.12) 


6 7 


B 21* G 12® B 22 * C 22 “ 1 (3.7,13) 

and 

B 21* C ll@ B 22 ,C 21 = 0 

The matrix B 2 may be chosen arbitrarily, with the condition 
that the matrix B’ is full-row ranked . Formats of matrix B, 
always suggests that the LIN vector of 8 2 may be obtained 
as [®2± r ^22^ ky choosing B 2 ^ as the null matrix and B 22 as 
the identity matrix of appropriate order, i.e. 


0,.Q© 1. 1 = I is satisfied 

and hence the resulting modal matrix takes the form M = [j], 
where Q is of the order n/2, n being the order of the Heed- 

tVjtX \ 

Muller matrix R and Q also takes the same form as the/matrix 
R, and I is of the order n/2. The identity that 

Rank + Nullity = Order is also verified as a result 
of this* 


The matrix Q in the modal matrix M, may also be 

defined as the matrix Q, the characteristic matrix, with 

Y«X, and is of order n/2 or the same may be obtained as the 

characteristic matrix, with y=l, for matrix K’ , obtained as 

(m-1) 

a direct-product of the core matrix times and I is the 

identity matrix of order n/2. Alternatively, Q may be written 
[ K * i]. The span of the column vectors of the modal-matrix 



I 
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as the basis v-ctors, the complete set of eigen-vectors of 
the RM matrix may be found. 


Illustrative Example 3.7.3: 

Given that the no. of variables in the switching 
function, m=2. 


The for the expression 


f(x 1 ,x 0 ) = c 0 to© c 1 t 1 © c 2 \© c 3 t 3 
where t's are the minterms and c ± (0,l\, indicates whether 

X 

the function value is zero or one for the related input 
and as the mintermsare disjoint, the inclusive OR, 

operation is replaced by the exclusive OR operation. For 
this example let the CjS be 01 01. The other sets of c ± s 
of all possible 2 m values that this may take are the 
sixteen combinations of Is and Os, from 0000 through 1111. 


R = 


The transform matrix, R for this example is 



”l0 00 ” 


"lO 00 

11 00 


00 10 

10 10 

OX 

11 00 

11 11 

I 

00 11 _ 


and the eigenvalues are 

1 i.e. y-l, 
1 


% 


= 4 
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To find the modal matrix, 
of characteristic matrix, 



we like to find the null space 
which m the simplified form is 


[° °] and I = fl 0, 




and 


M = 


0 0 
1 0 
1 0 
0 1 


and the set of all eigen-vectors 

r oooo 
0 10 1 
0 10 1 
0 0 11 


H = 


The characteristic matrix of the transform matrix R, 

0 0 0 0 


is (R-yl) = 


10 0 0 
10 0 0 
1110 


and then 


find the null-space the matrix B so obtained 


i.e* B = [J- ? ? £] and obtain the set of x. such that 
1 1 x U 


Bx = 0 
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and B * 


0 0 o' 

► 

X 1 ^ 

r 

0 

1 1 0_ 


x 2 


0 


X 3 


0 


X 4j 


0 


B. 


B, 





1000 

may 

be obtained 

as 

1110 

0010 




0001 

o 

H* 

O 

= [c L c 2 ] , 

we may write 


"lOOO* 

1110 


C 11 c 12 


10 0 0 

0 10 0 

0 0 10 

0 0 0 1 

m 

C 21 C 22 


0 0 10 

0 0 0 1 






'12 L 11 


ri°i + c r 00 i = r 00 i 

Lt t J * c 22 l 10 j l 00 J 


and c 12 [qq] + c 22 [ Q1 ] ~ [ Q1 ] 


c * [ 10 ] 
c 22 L 01 J 


and hence c l2 [^J + [q^J = 0 

rlO") rlOi rOO-j „ flOn rOOi = rOOi 
C 12 = ‘•11-' L 01 J L 10 J L 11 J L 10 J L 10 J 

(as the matrix [^] is involutary 


C 2 ~ 


00 

10 

10 

01 
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Illustrative example 3.7.4: 

the 

The number of variable in j_ positive canonic F»M 
expansion is given as m=5, 

a) The order of RM matrix is n = 2 m = 32 

b) The modal matrix may be easily constructed as 


OOOO 0000 0000 0000 
1000 0000 0000 0000 
1000 0000 0000 0000 
1110 0000 0000 0000 


M=[?]= 


1000 0000 0000 0000 
1100 1000 0000 0000 
1010 1000 0000 0000 
1111 1110 0000 0000 

1000 0000 0000 0000 
1100 0000 &000 0000 
ioio oooo aooo oooo 
1111 0000 1110 oooo 

1000 1000 1000 OOOOi 
1100 1100 1100 1000 
1010 1010 1010 1000 
1111 1111 ml 1110 

1000 oooo oooo oooo 
0100 oooo oooo oooo 
0010 oooo oooo oooo 
0001 oooo oooo oooo 

0000 1000 oooo oooo 
0000 0100 oooo oooo 
OOOO 0010 oooo oooo 
QQQO 0001 oooo oooo 

0000 0000 1000 oooo 
0000 OOOO 0100 oooo 
0000 0000 0010 oooo 
oooo 0000 0001 oooo 

oooo oooo 0000 1000 
oooo oooo OOOO 0100 

OOOO OOOO OOOO 0010 
OOOO OOOO OOOO 0001 



CHAPTER 4 


GENERATION OF FEED-MJllER EXPANSION COEFFICIENTS 

Developing algorithm , which synthesise any arbitray 
switching function, using a minimum number of logic circuit 
has been a perpetual problem in switching theory-* Boolean 
realisation use ANE/OR/NANE and NOR gates, whereas these 
functions may also be implemented using Exclusive OR gates* 

The advantages of designing with Exclusive OR gates, were 
outlined in Sec. 1.1* Applications using Ex-OR gate reali- 
sation has suffered because there are not ma#y algorithms 
that may be easily applied to find the minimum solution, as 
are available in the case of Boole a* vertex realisations. 

In the exclusive-OR realisation for the minimum design, it 
is desirable to minimise the number of terms, in exactly 
the same manner that it is desirable to minimise the number 
of prime complicants (AND ) terms in a conventional BooleiA 

realisation* 

Bigen structure of any matrix characterises the 
matrix and therefore it is possible to associate the eigen- 
structure of the Bee 4-Mu iler matrices, obtained in tbs previous 
ohapter, in the generation of the positive canonic- and non- 
positive eanonie Reed-Mullsr expansion coefficients . Thus 
the eigenstruoture of the expansions 1* able to relate the 
minted coefficients (obtained freothe truth table# of the 
Boolean function) to positive canonic and non-positive 
canonic RM expans ion coefficients. 
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The algorithms developed h<=re are for the two different 

RM matrices* In the algorithm, for the positive canonic RM 

h 

expansion, eigenvectors parta^ns the complete set of minterm 
coefficients vector space C, into cosets and "by associating 
the coset leader with an eigen vector, the positive canonic 
RM expansion coefficients may he obtained* The other algorithm, 
makes use of the eigenvalue of the non-positive canonic RM 
matrices and is able to associate the minterm coefficients to 
its corresponding non-positive canonic RM expansion coeffi- 
cient for a given pblarltjrf unction of the Boolean function* 

Both the algorithms may he programmed, when the matrices 
cannot be handeld i.e* for Boolean functions with large number 
of variables, and may he employed profitably in the minimisa- 
tion problems . 

4-1 PURPOSE OF GENERATION OP THE NON-POSITIVE CANONIC RM 
EXPANSION COEFFICIENTS 

The generalised Reed-Mu lie r (GRh) expansion may be 
expressed as (1.7) 

f<V *2' —• * V = a o + a 1 *1 + a 2 *2 + 

+ a 2 m -i ^ ^ 

where each 5^ input variable may appear in positive or negative 
form but not in both forms. 
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The number of zero-valued a^ coefficients in the 
GEM expansions for any given function f(x) and hence the number 
of remaining AMD and Exclusive terms m the resulting reali- 
sation, will vary depending upon the chosen polarity of each 
x^ input variable* The extreme polarities of positive and 
negative do not contain more serovalued terms than some other 
non-positive canonic RM expansions. Hence the criterion -for 
minimisation is the selection of that non— positive canonic 
expansions which has the maximum number of zero (minimum 
number of non-zero) a i coefficients in the expansion. A com- 
parison of all possible non-positive csncnic RM expansions are 
involved in order to determine a minimum non-positive canonic 
RM expansion [21], [22], [ 23 ], [27]. It is difficult to 
generate these 2 m positive and non-positive canonic RM ex- 
pansions and these are not intutive either. Constraints 
on the variable being complemented may give rise to minimum 
solution, A method has been developed [7], where geometric 
plotting and grouping have been outlined and thus after manu- 
pulating the coefficient values, a minimised non-positive 
canonic RM expansion may be achieved without involving 
exhasutive search procedure. A computerised in place algo- 
rithm has also been devised [10]. The algorithm enables 
sequentially generate all non-positive canonic RM expansion 
coefficients, and the algorithm being fast, an exhaustive 
seareh for the minimum is feasible. 
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The table below clearly shows the advantage of the 
non-positive canonic RM expansions. 


TAB I® 4.1.1 


EEQUlRBi^NTS FOR POSITIVE AND EON-POSITIVE 
CANONIC RM EXPANSION IMPLEMENTATION 


s. 

Number 

Minterm* 

Positive 

Non 

-Positive 

— 

Polarity* 

No. 

of vari- 

Coeffs . 

| Canonic 

RM 

Canonic 

RM 

Pattern 


ables in 
the 

f{x 1 x n 

j)B2tpans ion 

Expans ion 

of the in- 
put vari- 


function 


X.0R 

AND • 

SOT 

TOR 

AND 

NOT 

able 

f 

Mm 

L~4 

.2228 

6 

5 

0 

2 

2 

2 

6 

2 

4 

21845 

2 

1 

0 

1 

0 

1 

1 



(101010...) 








3 

4 

65280 

1 

0 

0 

1 

0 

0 

0 



_LQ -*^001 ...11^| 







5 

4 

4 

43690 



1 

0 

0 

1 

0 

0 

0 



(0101 ...) 








_5__ 

4 

255 

2 

1 

0 

1 

0 

1 

8 

6 

5 

1 OQOO ! 

16 

14 

0 

6 

m 

MB 

\ IT 



i 




5 

B9 


| 23 

_7_. 

5 

20000 

ifL 



_12 

_2_. 

__6_ 

H 

■ 

! 16 

8 

5 

1421655765 ^ 

_2_ 

i 1 

0 


0 

1 

i 1 

mmm m 

■9 

5 

227216640 

.16.. 

um 


7 

h _6__ 

_2 

14 

Eq__ 

5 

1671935 

_2_, 

1 

0 j 

n 

0 

_1__ 

8 

m 

6 

1 1 0000000000 

22_. 



m 

BH 


54 

IfPJ 

6 

.21241814220 

24- 

. 21 .. 

0 

12 

12 

i_4 

60 


6 

.22121694128 

2§- 

,22__ 

0 

14. 

_14__ 


61 

SB 

6 

.1224507890 

20_ 

18 

0 

10 

10 

1 

22 

U_ 

6 

.1020204920 

2§_ 

.24__, 

_0__ 

11_ 

_il__ 

_2 

42 — _ 


* Decimal value of the binaryrepresentat ion of the minterm co- 
efficients and the polarity function. When the pattern of rain- 
term coefficients is not indicated, then it is a random pattern. 
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Minimisation criterion, depending on the realisa- 
tion may "be applied ; one is to minimise th° number of AND 
gates, another one is to minimise the overall number of inputs 
to the Ex— OR gates or to minimise the overall number of 
inputs. However, the choice of criterion being the economical 
imp 1'- mentation of the situation to which the design is to b c 
applied. The penalty to be paid for the non-positive canonic 
RM expansion implementation is th^ sub-optimum number of 
product terms, while the number of input variables or the 
number of inverters required is gen -rally less than for the 
minimum form [$0]. 

4.2 GENERATION OP POSITIVE CANONIC RM EXPANSION COEFFICIENTS 
USING EIGENVALUES 

The vector*' space of all combinations of possible 
values that c ± s may take get partioned by the eigenvectors of 
the transform into subgroups, and the RM coefficients for a 
particular subgroup are governed by a particular eigenvector* 
Thus it is observed that this vector space is partioned into 
cosets and the elements of the subgroup or the cos^t is in- 
variant i.e. the RM expansion coefficient is found from within 
the coset, governed by an eigenvector. 

This if h t , h 2 , .... h* are the eigenvectors of the 
transform, then cosets may be formed as shown (all operations 
are modulo 2 addition) 



77 


i 2 3 

k-j S-j (i)^ g i 

e 2 0 h 1 g 2© h 2 S 2© ll 3 "•• 



g k© h 1 g k © ^2 g k© h 3 •*•• s k© h n 

and in this all g 1 (+)h^ s become to coset leader. While 
forming the cosets, instead of g^Qh^ being any element that 
has not appeared earlier, it is the n<=xt lowest element after 
h.j , that has not appeared earlier. Thus the coset leader is 
the minimum vector m the coset. In the algorithm bit-by-bit 
modulo 2 addition of all eigenvectors to the c^s gives the 
entire coset and by chocaing the minimum vector the coset* 
leader is thus obtained. 


Now after identifying tht coset leader, the eigen- 
vector that relates the c^s to the RM expansion coefficient 
may he found by multiplying the coset leader with the charac- 
teristic matrix of the transform. Having found the eigen- 
vector that relates the c^s to the Rh expansion coefficient, 
a mere bit-by-bit addition of o to the eigenvector gives the 
RM expansion coeff icaents . 



i 

Qi * 

W 

PI V: 




OJ 

* 


PH 


The algorithm is schematically depicted above. 
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Illustrative Example 4,2*1 

The cosets for m = 2 could he written as shown below 
(the decimal value of c^s are given for ease of notation)* 


Coset 0 

1 .Leader”"* 

3 

2 

i 

^ 


0 

6 

8 

14 

Eigenvectors 

4- 

1 

1 

7 

9 

15 

14 

2 

2 

! 

4 

10 

12 

8 

3 

3 

5 

1 1 

13 

1 

6 


Thus for the given c^ 0101, which is notationaly 10, 
(Ten) by car - ying out the bit-by-bit operation on all eigen 
vectors we get 10, 12, 2 and 4, the minimum of which is 0100, 
which is 2, the coset leader also* Now carrying out a vector 
matrix multiplication, i.e. the coset leader and the charact- 
eristic matrix, 


”0“ 


“o 

0 

0 

0 ! 


’« ! 
t 

0 


1 

0 

0 

0 


1 ! 

0 

j 

1 

0 

0 

0 


0 

1 

1 i 


1 

1 

1 

0 ! 


0 

i_ _ 


we obtain the RM expansion coefficient determining eigen 
vector, which is 0001, and bit-by-bit addition of this to the 
c^s gives the Rl-'i expansion coefficient. 
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i.e. 


r 

I 0 

I 0 


! 




+ 



1 ; 

: which, is so. 


Illustrative Example 4,2.2 

The coset formation for a 3 variable function, is shown 
in the table. The coset formation contains all the 2 n 
possible combination which the set of c ± s , the minterm coeffi- 
cients could take. Hence once the table is formulated and 
the appropriate coset leader associated with th° corresponding 
eigenvector, using the algorithm} the table is available as 
a ready reference or as a mere look-up table. The effort 
therefore is ouly to do a bit-by-bit addition, modulo 2, of 
the minterm coefficients with the corresponding eigenvector 
to give the Reed-Muller expansion coefficients. However, 
in this table the positive canonic RM expansion coefficients, 
corresponding to each minterm coefficient is also given as 
the denominator to it i.e, after the vector obtained as a 
result of bit-by-bit addition of the minterm coefficient 
with the positive RM expansion coefficient determining eigen- 
vector found from the table. 



TABLE 4.2.1 

G OSE'T XBADE R. BIG-EiWECTOR LOOK-UP TAB LE 



The above table p.lso gives the positive canonic RM expansion coefficisnts 
The denominator of each of the minterm coefficient entry is the posijive 
canonic RM expansion coefficient corresponding to it. 
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Illustrative Example 4 , 2*3 

The number of variables in the positive canonic RM 
expansion is given as m =* 4, 

(a) The order of the RM matrix is n = 2 m = 16 

(b) The modal matrix may be obtained as 


R = [«] = 


00000000 
1 0000000 
1 0000000 [ 
1 1 1 00000 
1 0000000 
1 1 001 000 

10101000 i 

11111110 I 

1 0000000 . 
01 000000 I 
001 00000 ! 
00010000 
00001 000 
00000100 
oooiooi 0 
00000001 


The above linearly independent eigenvectors as the 
basis, the span of these vectors may be found to be consis- 
ting of the 256 eigenvectors of the RM matrix. There are 2 n 

1 6 

possible combinations of the minterm coefficients i,e. 2 =6553^ 

So instead of handling this many combinations, we may store only 



the 256 eigenvectors of the RM matrix and hy "bit-by-bit addi- 
tion of the given minterm coefficient to these 256 eigen- 
vectors and picking the minimum vector, the coset leader is 
obtained* A 

The process of associating the expansion coefficient 
determining eigenvector to the coset leader could be available 
as a look— up table, as this is independent of any operations 
with minterm coefficient and therefore once the cosct leader 
is obtained, then by identifying the expansion coefficient 
determining eigenvector, bit-by-bit addition of the minterm 
coefficients gives the corresponding RM expansion coefficients, 

This algor ithm is very efficient, provided the two 
look-up tables are available before hand viz* the set of all 
eigenvectors and the coset-leader to eigenvector table, then 
with only two sets of bit-by-bit additions, the RM expansion 
coefficient may be obtained. In general, 2 n//2 + 1 bit-by-bit 
addition and the process of picking the coset leader is ade- 
quate and the coset leaders could also be stored in memory 
for comparision and thus speed-up the operations. The number 
of coset leader being the total number of eigenvectors itself. 

4*3 GENERATION 0? NON-POSITIVE CANONIC RM EXPANSION COEFFICI-*: 

ENTS, USING EIGENVALUES 

To enable relatathe o^s to the positive and non-posi- 
tive canonic RM expansion an algorithm has been developed, 
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which is schematically shown below 





1 

Jt 

<1 

— jy~ 

L 

Characteristic | 
■matrix i- 

i 

¥ y s 


{ 

i 

i 


Goeff, 
jM deteraun 
mg . 

n.a.c:kpr__ 


-4 


*j~\ — — — 


GRM 

coeff ic lent 


This algorithm may be seea as a version of the one 
presented in the earlier section and this being applicable 
to the more general case may be applied to the ’positive 1 
polarity RM expansion coefficients also. 

This algorithm may be written as 

[| y T (R + y I)] -[X] + [X] = RX 

There are three kind of eigenvalues that are obtainable from 
the RM matrices viz. y = 1 , y = a and y = 1 + aj also it was 
mentioned in section %2 that whenever y = a occurs then 
y = 1 + a also occurs. The total number of eigenvalues is 
2 m and thus the number of eigenvalues that appear for y = 1 , 
if they occur, is always even. Bft must be remembered that 
all arithmetic is modulo 2 and therefore the summation of 
y(R-yl) for all ys , as far as y = 1 is concerned may be omitted 
and evaluating for y = a and y = 1 + a, we obtain. 


iHS = [ a(R + al) + (r + « )(B + (1 + a) I)][x] + [x] 
- [aR + a 2 I + (1 + a) R + (1 + a)^ l][x] + [x] 
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- C 0 ®- + ! + «I + E+ c®. + al] [x] + [x] 

- [i + R] [x] + [x] 

= X + RX + X = RX = RHS QED. 

It is observed that "weighting the characteristic 
matrix "'toy the root results in the characteristics ^matrix, 
which is what has been used in the previous algorithm. 

In this algorithm instead of obtaining the eigenvectors that 
determines the RM expansion coefficients, some arbitrary 
vector is obtained, which determines the GRM expansion coe- 
fficients. 


Illustrative Example 4 * 3*1 


For a 3 variable Boole an, function of polarity x^ x 2 x^ 
i„e, m = 3, polarity = 101 


R 



x 


[ 01 ] 
L 11 J 


X 


C 10 ] 

L 11 J 


fooi ooooo 
001 1 0000 
1 01 ooooo 
1111 0000 
001 0001 0 
001 1 001 1 
1 01 01 01 0 
11111111 


1 0000000 
001 ooooo 
00001 000 
0000001 0 
1 1 000000 
001 10000 
00001 i$>0 
0000001 1 


01 000000 

0001 0000 

00000100 

00000001 

1 1 000000 

001 1 0000 

00001 1 00 

0000001 1 
«•» ■*“ 


1 0000000 
001 ooooo 
00001 000 
0000001 0 
11000000 
001 1 0000 
00001 1 00 
- 0000001 1 


for the given set of o* say 10011001, we may obtain the GRM 


expansion coefficients as 


a = R " c 
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Alternate ly t we know that there is only one variable being 
complemented i*e t = 1 and we may obtain the number of roots 
in the ground field and extension field either as direct 
product of the individual core matrices such as 

[j] » t, *„] * t, 1 ] 

i.e. r = * ; »v = 4 

v = 1 + “ i m 1+a = 4 

or using the formula we have 

m = 3 , t = 1 therefore there are no roots in the 

ground field and all roots are in the extent ion field and 

there are d = 2 m of them i.e. 8, 4 of them are y = <x and the 

2 

other 4 are y = 1+a (or a ) 

The characteristic equation Is 
2 

(y + y + 1 ) =0 

i.e, (y 2 + y + I) 4 = 0 

Now using the algorithm the G-E.M expansion coefficients 
can be found as follows : 

a = [£ y(R $ (£> Cc] 
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\ 


a 

0 

1 

0 

0 

0 

0 

0 


1? 

0 

1 

0 

0 

0 0 

0 



1 ~ 

C 1 

• c i 

0 

a 

1 

1 

0 

0 

0 

0 


0 

a 2 

1 

1 

0 

0 0 

0 




i 

G 

















2 


2 

1 

0 

a 

2 0 

0 

0 

0 

0 


1 

0 

a 

0 

0 

0 0 

0 



Cry 




















3 


3 

1 

1 

i 

a 2 

0 

0 

0 

0 


1 

1 

1 

a 

0 

0 0 

0 



c „ 


c 









U lea 










4 


4 

0 

0 

i 

0 

a 

0 

1 

0 


0 

0 

1 

0 

a 2 

; o t 

0 

f 

i 1 

1 

t 

1 

w 

c 5 

0 

0 

i 

1 

0 

a 

1 

1 


0 

0 

1 

1 

0 

a 2 1 

1 


i 

I 

1 

°6 ! 


c 6 

1 

0 

i 

0 

1 

0 

o?0c 

I 

1 

0 

1 

0 

1 

0 a 

0 

1 

i 

* 7 ; 

i 

i 

i 

i 

C 7 

1 

1 

i 

1 

1 

1 

1 

a 


1 

1 

1 

1 

1 

1 1 

c 



i 

_ C 8J 

i 

i 

C 8 ! 








— M*“. 







— 

J 


1 


and simplifying over the irreducible polynomial, 


~~1 01 00000 ~ 


V 


°1 ‘ 

01 1 1 0000 


c 2 


°2 

1 0000000 


°3 


°3 

1 f 1 00000 


°4 

© 

C 4 

00101 010 


°5 


c 5 

001 1 01 1 1 


c s 1 


c 6 

1 0101000 

! 

c 7 

i 

c 7 

11111110 


c 8 

1 

t 

_° 8 _. 


“e^ 

c 2 ec ^ 4 


c t ac 2 ^ 5 

c^^Cpjecy 

°^ c 4 ^^ c 7^ c 8 

c 1 @c^c 5 

c 2 ^ c ^ 5 ^° g>^ c 7 



■° 1 1 


C 2 


°3 

0 

C 4 


c 5 


°6 


C 7 

) 

C 8 



winch for the given set of c s viz. 10011001 , is 
a T = [01100000] 

and thus the required GRM expansion coefficients have bpen 
found for the polarity 101. 

4.4 GENERATION OE THE Rh EXPANSION COEEBICIENIS BASED ON GOOD'S 
E ACTOR IS AT ION OE THE RM MATRICES 

An *m» variable Boolean function could have any 
polarity for specified optimum solution. As already mentioned 
•in-place 1 computer algorithm [10] is available for this. 

The program given at Appendix & t©o, though not 'in-place 1 , 
gives the same results, with the flexibility that even for any 
random polarity the non-positive canonic RM expansion coeffi- 
cients may be had in a single run. The Reed-Mu Her matrix 
is built-up as Kronecker product and then factorised as two 
kinds of matrices using Good’s technique, which when sequen- 
tially processed gives the desired pos it ive /non— posit ive canonic 
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RM expansion coefficients. As the representation in the 
program for the matrix is the existence or otherwise of an 
»1 » element m the matrix, the computation is also fast i.e. 
only the location of the elements in the matrices determines 
the output vector. The Reod-Muller mebrix is thus factorial 
as 




with R. taking one of the two formats depending on wheth-r 
the variable is true or complemented. The program in its 
present form scans all the possible 2 polarity functions in 
a binary sequence and from this the minimum number of AND gates 
or inputs to Exclusive OR gates or with any oth*r specifications, 
the desired optimum polarity may be had. Ihe limitation of 
this program is the capability of the computer/ handle matrices. 
The program is given at Appendix £> • 


The minterm coefficients for a 5 variable Boolean 

-t 

function is given as the input. These minted coefficients 
are the same used in [10] for the example there. The output 
may he had in two forms, one as a summary of the number of 
gates required for each polarity function which is reproduced 


as a table her®* 
gives the details 


the other besides giving the ebove data 
of RM expansion coeff icisnts , a x s for each 



case. The program in its present form does not give the 
exact number of Ex— OR gates required, though the same may be 
construed and neither does it indicate the type of the AND 
gates required i.e* how many input AND gates. Howev^r^ an 
analysis of a^s, the RM expansion coefficients would divulge 
the number of inputs to each Mb gate. 

TABLE 4-4.1 


RM EXPANSION COEFFICIENTS FOR DIFFERENT POLARITY OF THE 

VARIABLES OF A BOOLEAN FUNCTION 


B-J 
Super- 
s cript 

Number of 
input to 
ExOft gates 

Number of 
AND gates 

Number of | 

Inv - rtors 1 
(NOT gates) 

Re marks 

0 

16 

14 

0 


1 

14 

13 

1 


2 

12 

11 

1 


3 

14 

1 2 

2 


4 

9 

9 

1 


5 

1 1 

1 1 

2 


6 

8 

8 

2 


7 

13 

12 

3 


8 

17 

14 

1 


9 

16 

14 

2 


1 0 

12 

10 

2 


1 1 

16 

13 

3 


1 2 

9 

8 

2 


13 

10 

10 

3 


14 

7 

6 

3 


15 

11 

10 

4 


16 

20 

17 

1 


17 

18 

15 

2 


18 

18 

15 

2 i 


1 9 

16 

14 

3 ! 


20 

13 

12 

2 


21 

17 

15 

3 


22 

14 

13 

3 


23 

21 

17 

4 


24 

1 6 i 

14 

2 


25 

16 

13 

3 


26 

14 

12 

3 


27 

14 

*2 

4 


28 

10 

9 

3 


29 

12 

12 

4 


30 ; 

10 | 

9 

4 


31 

14 

12 

5 
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4.5 GE HERAT I OlsT OF MIHTERI v i CCEFFICIEHIS FROiri A G I'VEH Rk 
EXP ANSI OH CCEFFIOISHTS , FOR TEE GIVER POLARITY 

The linear transformation from the mint^rm coeffi- 
cients to the RM expansion coefficients was represented by 
the RDu matrices . The expression being 

a = Rc (4.5.1) 

Therefore from the given set of non-positive canonic Rk 
expansion coefficients, the minb^rm coefficients may be 
obtained as 

c = R H a (4*5.2) 


The RM matrices were obtained as direct product of core 
matrices, and structure of the inverse of the RK matrices is 
also obtained as direct product of the inverses of the core 
matrices. The inverse of the core matrix is the core matrix 
itself, the core matrix being involatary. The mv'-rse of the 
core matrix, for a complemented variable is 


r° I]" 1 = [} \] (4.5.3) 

L 1 1 J L 1 0 J 

and therefore the inverse of the RM matrices are oh tain da s , 
direct-product of these inverses comb me d as per the polarity 
function of the input variables * Thus for a polarity 
function 101001, the R ^ is 
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Illustrative Example 4,5.1 

J?or tho 3 variable Boolean function consider?! 
earlier, with polarity function HO, the non-positive canonic 
RM expansion coefficients were obtained as 

[a] T = [01 100000] 

Iho matrix R~ 1 may be obtained as 

R -1 = [ } °] a t] °] & [J i ] 

Therefore the nunterm coefficients may be obtained as 


1 1 000000 


0 ! 

1 0000000 


1 ! 

1111 0000 


1 1 

1 oi ooooo ! 


0 ! 

1 1 001 1 00 


0 * 

1 

1 0001 000 


0 i 

11111111 


0 1 

\ 

10101010 


0 1 
_ i 


which i$ [c] T = [10011001] 

These were the minterm coefficients originally assumed 
in the example and therefore r" 1 may be obtained as mentioned 
in this section. 



CHAPTER 5 


CONCLUSIONS AND SUGGESTIONS FOR FURTHER RESEARCH 
5 , 1 SUMM JIY OF THE THESIS: 

Minimising the number of logical circuit elements 
has always been a problem. Various techniques such as 
Karnaugh map and Ouine McClusby’s technique, and the 
recent 'b mapping’ help achieve a limited aim. The 
search was for more economical realisation and therefore 
designing only with exclusive-OF* and ,.ND cates were 
adopted, ( in addition, mvertgsME too would be required). 

This gave rise to the problem of minimisation of the Lx-OR gates 
gates and the parameters on which this minimisation has 
been defined is the polarity function of the variables of 
the boolean function. So the boolean function, via its 
minterm coefficients could be related to which came 'to be 
known as the hM expansion coefficients. The transformation 
was accordingly termed as the Rivt matrix. 

The RM matrices were obtained on direct-product of 
core matrix and thus the eig^n-values could be easily 
obtained as direct— product • Ihe problem of multiple roots 
of the characteristic equation, especially, m the case o± 
the RM matrix for the positive canonic RM expansion was 
simplified and the eigen-values of the RM matrices for the 
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non— posi tive canonic RM expansions were formulated, aqain 
as a function of the polarity function of the input variables. 
As the RM matrices could be obtained as direct-product of 
core-matrices, they could also be factorised by Good's 
technique. The main advantage of this factorisation has 
been that the few non-redundant entries in a matrix could 
be made fewer, resulting m fewer number of operations for the 
transformation implementation. This aspect has been exploited 
utmost for computer implementation of the generation of 
the RM expansion coefficients. 

The purpose of generation of RM expansion coefficients 
itself was considered essential for minimisation of the 
Ex-OR designs. The eigen-structure of the RM matrices were 
used for the generation of the RM expansion coefficients. 
Good's technique of factorisation was also put to use. 

The main dark area in the use of Ex-Oh design has 
been that not many usable algorithms are available. This 
thesis is an effort to contribute towards understanding the 
problem and try and give an easier approach to it. Table 
4* 1 ,1^/Chapter 4 would evince tho circuit designers into 
thinking of designing chips, thaugh may be exclusively 
dedicated, but extremely small in size. However, the 
economical production of the 'chips' would be the deciding 


factor 
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5.2 SUGGESTION FOR FUTHLh RESEARCH: 

The eigenvectors of the Ri a inacric.es, for the positive 
canonic RM expansion, could be used to form lookup rabies 
and thus by a mere bit-by-bit addition, of the mint arm 
coefficients with the expansion coefficient determining 
eigenvector, the expansion coefficient could be obtained. 
Similar lookup tables could be formed for the RM matrices, 
for the non-positive canonic RM expansions, so that the 
non-positive RM expansion coefficients too could be obtained 
m a simple manner. 

Practical designing of the circuits, using cxclusive-OR 
gates, AND gates and invertors could be profitably 
worked out. The algorithm developed in this thesis or 
ClO ] t could be put to use. 

Lastly, every digital information processing essentially 
being represented by a set n iunction of m vaxiables on the. 
GF(2), such a set of functions could be reduced to one 
polynomial of one variable over the extension field, GF(2 N ) . 
Such polynomials have remarkable properties, which may 
be used effectively in designing of "itching circuits [28]. 


<1 



REFERENCES 


1* Murago S., ’Logic Design and Switching Theory’ 

John Wiley and Co., 1979. ' 

2. Chinal J», 'Design Method for Digital Systems', 
Springer-Verlog, Berlin, 1967. 

3, Reed IS, ' A class of Multiple— error Correcting Code and 
the Decoding Scheme', ire Trans. 1954, IT-4, pp. 38-49. 

4* Muller DE, 'Applications of Boolean Algebra to Switching 
Circuit Design and to Error Detection', IRE Trans. 

1954, EC-3, pp. 6-12. 


5. Ashenhenhurst RL, ’The Decomposition of Switching 
Functions', Proc. Intelnl Symp. on the Theory of 
Switching, Harvard University, 1959, vol. 29, pp. 74-116. 

6. Roth and Karp, 'Minimisation over Boolean Graph', IBM 
Journal Ros. and Development 1962, vol. 6, pp. 227-238. 


7* Wu, Chen and Hurst, 'Mapping of Reed-Muller Coefficient 
and the Minimisation of Exclusive-OR Switching Function' , 
IEE Proc. 1982, vol. 129, Pt. E, Mo. 1, pp. 15-20. 

8. Reddy SM, 'Easily Testable Realisation for Logic Functions', 
lEElTrans * 1972, C 22, pp. 1183-118 


9. Birkhoff and MacLano, 'A Survey of Modern Algebra', 
MacMillan and Co., 1962. 


10. Besslich P.D , , 'Efficient Computer Method for Exclusive-OR 
Logic Design', IuE Proc. 1983, vol. 130, Pt. E, No. 6 

pp. 303-306. 

11. Antan Glaser, 'History of Binary and other Non-decimal 
Numeration', published oy the author Southampton, Penn. 
USA. 


12. Bellman R., 'Introduction to Matrix Analysis' McGraw Hill, 
New York, I960. 


13. 


DeRusso, Roy and Close, 
John Wiley and Sons Inc 


'State Variables for Engineers 
1967, PP. 255-262. 


14. Kurosh A., 'Higher Algebra ', (English Translation) 
Mir Publishers, Moscow. 


f 



97 


15. 

16 . 

17. 

18. 


19. 

20 . 

21 . 


22 . 

23. 


24. 

25. 

26. 


27. 

28. 


Lanchanster P., 'Theory of Matrices’, Academic Press, 

Now York, 1969. 

Pease jviC, 'Methods of Matrix rt lgcbra' academic Press, 

New York, 1965. 

Peterson WW, 'Error Correcting Codes', a. IT Press, 1961. 

Samuel u. Lee (edited by MNE DC), 'Computer Science and 
Multilevel Logic Theory and Applications North-Holland 
Publishing Co., 1977. 

Searle SK, 'Matrix algebra useful for Statistics', 

John V/iley and Co. 1982. 


Williams, CL, 'Boolean nlgebra with Computer Applications', 
MacGraw Hill, 1970. 


Mukhopadhyaya and Schmitz, 'Minimization of Exclusive-Oi- and 
Logical Equivalence Switching Circuits', ISLE Trans. 1970, 
C-19, pp. 132-140. 


Papakonstantmov G., 'Minimisation of modulo 2 sums-of- 
products', IELE Trans. 1979, C 28, pp. 163-167. 


Kodandapam and Setlur , 'A note on Minimal Reed-Muller 
Canonic Forms of Switching Functions', IEEE Trans. 1977, 
C 26, pp. 310-313. 


Hurst SL, 'Logical Processing of Digital Signals', Edward 
Arnold, London and Crane Ruesak, New York, 1978. 


Andrews and Caspar! , *a General Technique for Spectral 
Analysis', IEEE Trans. 1970, C-19, pp. 16-25. 


Good IJ, ' The Interaction Algorithm and Practical Fourier 
Analysis', J. Royal Statistical Society, London, B— 20, 
p. 361,. 1958. 


Saluja and 
Expansions ' 


Qng 'Minimisation of Reed-Muller Canonic 
, IEEE Trans. 1979, C-28, pp. 535-537. 


Takahashi I., 'Switching Functions Constructed by 
Extension Fields', Information and Control, 1981, vol 
48, pp. 95-108. 





